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TWO DIMENSIONAL THEORY OF STIFFENED 
PLATES 


by 


CHANG O’CHOU, Ph.D. 
(The University, Bristol) 


SUMMARY 


The two-dimensional stiffened plate problem is completely governed by the 
fundamental differential equation 


where k, and k, are two non-dimensional structural constants depending on the 
relative proportions of the reinforcing members. 


The complete solution is 
o=f, (x+ik,y)+f, (x—ik,y) +9, . (25) 


This solution renders all solutions for the isotropic plate available for the stiffened 
plate also. The choice between the several types of solution available is governed 
by the same considerations as for the isotropic or unstiffened plate. 


1. INTRODUCTION 


By the term “stiffened plate” is meant a plate of thickness t, reinforced by 
longitudinal members (stringers), each of section A,, at spacing b, and transverse 
members (ribs), each of section A, at spacing a, as shown in Fig. 1. Both longitudinal 
and transverse reinforcing members are continuously attached to the sheet, and 
are treated as of negligible stiffness in bending. Thus, along the line of attachment, 
the displacement is the same for the sheet and the reinforcing member. A rigorous 
analysis of such stiffened plates results in partial finite difference equations, the 
solution of which is not familiar at the present time. However, the general stress 
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distribution can be obtained by considering both sets of reinforcing members as 
uniformly distributed; in other words, as two hypothetical sheets with thickness 
t,= A,/b for stringers and thickness t,= A,/a for ribs. The hypothetical sheet differs 
from the actual sheet in that it can only transmit normal stress in the direction of 
the reinforcing members which it represents and it is discontinuous in the direction 
perpendicular to the reinforcing members. 


Throughout this paper the x-axis is in the direction of the stringers and the 
y-axis in the direction of the ribs, as in Fig. 1. 


Notation 


A,, A, area of section of ribs and stringers respectively 
spacing of ribs and stringers respectively 
rectangular co-ordinates (see Fig. 1) 
displacements in x and y directions respectively 
see Fig. 2 
Young’s modulus 


Poisson’s ratio 


fz, fy, q normal components of stress, and shear stress in plate 


d, 1 dimensions defining the boundaries of a rectangular plate 
(see Section 3) 


rx /I 
thickness of plate 
A,/a 
A,/b 
defined in equations (19), (20) and (21) 
defined by equations (14), (15) and (16) 
see equation (13) 
stress function 

k,, k, defined in equations (23) and (24) 

k,’, k.’ defined in equations (29) and (30) 

K,, K., K;, K, constants of integration 
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X 


stringer 


rib 


2. GENERAL ANALYSIS 


Let u and v denote the displacements in the x and y directions respectively. Then 
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Ou 0 
Ox oy 


f. in plate= (1) 


(2) 


Nay * a: 


f, in stringer sheet = E ; (4) 


ov 


f.. fy, q are respectively the normal components of stress and the shear stress in 
the plate. 


Thus the total loading per unit width of stiffened plate in direction x is 


(= 
fs t= 


<4 Ou 
Oy. ox 
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and in direction y is 


Et Ou dv ., Ov 


while the total shear loading per unit width is 


By adopting these definitions of f, and f,, the conditions for equilibrium become 


Of. _ 
| 
Ofy oq 
and 0. 


fe (10) 
fy= (11) 


Substituting (10), (11) and (12) into (6), (7) and (8), eliminating any two of x, 
u and v, and using @ to represent the third, it follows that 


A, 

bt 

at 


Equation (13) is the fundamental differential equation for the stiffened plate. 
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As usual, a stress function . may be assumed such that 
a 


(7) 


(8) 


(9) 


3) 


4) 


5) 
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For problems with given boundary loadings, it is more convenient to treat 9 as 


representing x. In this case u and v are given by 
(17) 
where (19) 
on 


The similarity between (17), (18) and the familiar stress-strain relations for 
isotropic plates is obvious. The three non-dimensional constants z,, 2, and 2, 
may be called the x-ratio, y-ratio and xy-ratio of the stiffened plate. 


3. GENERAL SOLUTION OF THE FUNDAMENTAL EQUATION (13) 


It can be proved that the discriminant Q* — 4PR of (13) is always positive. Thus 
(13) may be rewritten as 


(5 ki? dy 2 
with . .  . (23) 
and 1/k,2=Q/2- . . . 


The general solution of (22) is 
(x+ik,y)+f.(x—ik,y) + 9, (x+ik.y) + (x—ik.y) & 


For stiffened plates with rectangular boundaries equation (25) may be expanded 
into double series with suitable periods, i.e. 


v k 


0 cos pé sinh k,’py p=0 cosh pé cos k,’pn 


(26) 
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kK : ‘ 29) 


and the four boundaries are x= +/ and y= +d. It is understood that A, assumes 
different values for different combinations. 


4. A SIMPLE EXAMPLE 


Consider for example a rectangular stiffened plate under pure bending in its own 
plane. Let 


x=y’, 


which is a solution of (13). Then 


Terms containing K,, K,, K,, K, represent rigid body displacements and can 
therefore be ignored. 


Both u and wv are single-valued everywhere; hence (31) is a solution. From 
(1), (2), (4), (5) and (12) it follows that 


f. in plate= ‘ BF 

f, in plate= (R-l)y . (38) 

f, in stringer sheet=62,y . (39) 

f, in rib sheet= — 62,,y ‘ . (40) 
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Along the boundaries y= +d, both f, and q vanish. There is therefore no 
resultant force acting anywhere along the horizontal boundaries. 


Along the boundaries x= +/, shear vanishes, while normal stress f, is 
proportional to the distance from the x-axis, i.e. the distribution is linear. The 
resultant force acting on the vertical boundary is zero while the resultant moment is 


d 


M= 6y*tdy=4td*. . . (42) 


d 


The stiffened plate is thus in a state of pure bending under this couple. The sign 
convention is such that a clockwise moment is positive at the right hand end. 


5. ANOTHER EXAMPLE 


Another example is a semi-infinite plate under a concentrated force acting 
perpendicular to the boundary. Let 


x= (x + ik, y) log (x + ik, y) + (x — ik, y) log (x — ik, y) — 
(x +ik,y) log (x +ik.y) —(x—ik.y) log (x—ik,y) (43) 


which is a solution of (13). Then 


2k,?x 
2x 2x 
2k,"y 2k."y 
= x74 k,?y? T x?4 (46) 
Eu= —(k,?2,+ log (x? + k,*y?) + (47) 


Ev=[2(a,+k,?2,,)/k,] tan“ (k,y/x) (48) 


In general, there is a line of discontinuity of v along the y-axis. Consequently, 
(43) is not a solution for plates containing the y-axis. However, both u and wv are 
single-valued either in the right half plane or in the left half plane of the y-axis; 
therefore (43) is a solution for the semi-infinite plane or any portion of it subjected to 
the edge stresses dictated by the stress formulae. 


Take the right half plane as an example; all boundary stresses vanish except 
at the origin which is a singular point. To evaluate the total force acting at the 
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x 


Yo) 
Fig. 2. 


origin, the condition of equilibrium of the portion of the plate including the origin 
is applied, as shown in Fig. 2. 


Because of symmetry, F is given by 


F=-2t wee) dx—2t\ f dy 


Therefore, this example corresponds to the classic radial pressure case for the 
isotropic plate and shows how the stress distribution is distorted by the presence 
of reinforcing members. For plates stiffened in the x direction only, this problem is 
usually classified under the heading “load diffusion.” 
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THE WHIRLING OF A SPINNING TOP 
by 


J. MORRIS 


(Consultant, Structures and Mechanical Engineering Departments, Royal Aircraft 
Establishment) 


SUMMARY 


This paper has been compiled because of the applicability of the treatment of 
the stability of the motion of the common spinning top to the problem of the 
whirling of rotating shafts, carrying loads of appreciable moments of inertia. This 
problem is assuming renewed interest and importance, especially in the drives of 
contra-propeller systems and the more recent high speed prime movers. 


1. INTRODUCTION 


This investigation deals with an old problem partly from a new angle. The 
reason for the resuscitation of the treatment of the stability of a common spinning 
top is because the analysis for the elucidation of whirling phenomena in rotating 
shaft systems is analogous to the classical analysis of the disturbed motion of a top. 
Whirling has become once again a problem of current and urgent importance 
primarily in pursuance of the recent development of light weight prime movers whose 
rotors revolve at very high speeds and secondly in connection with the mechanism 
for operating contra-propeller systems. Apart from these considerations the present 
treatment would appear to merit attention on its own account, as the question of the 
whirling of a top in consequence of unbalance does not appear to have been probed 
previously. It will be shown that the ordinary symmetrical top may have only one 
whirling speed and this will be above the minimum speed for stability in the vertical 
position. The unsymmetrical top however, in certain circumstances, may have two 
whirling speeds between which the unforced disturbed motion is exponentially 
unstable; and in general the conditions for its stability are somewhat more compli- 
cated than for the symmetrical top. 


Lanchester in Appendix VII of his “ Aerodonetics ” remarks 
“Tt is usual to treat the gyroscope and problems involving gyroscopic actions 


exclusively by mathematical methods in which the physical reality of the problem 
is partially obscured. One result of this is that, in spite of the unimpeachable nature 
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of the said methods, mistakes are not infrequently made, even the direction of the 
gyroscopic torque is sometimes misstated.” 


Discussing the gyroscope precessing in a horizontal plane, Lanchester gives the 
correct result for the gyroscopic torque, namely *=/m where /, » are respecti:ely 
the moment of inertia and the angular velocity about the principal axis, and ‘: is 
the angular velocity of precession. But when he proceeds to deal with the case 
when the axis of the gyroscope precesses at an angle z to the vertical, Lanchester 
derives the incorrect result namely, 7=/w{2 sin z, thus commiting a curious error 
in his physical reasoning concerning a basic proposition in gyro-statics. 


2. KINEMATICS OF THE PROBLEM 


The problem considered is that of a top spinning on a perfectly rough 
horizontal plane. 


Referring to Fig. 1, let O be the point of contact with the ground of the peg 
of the top, of which OC is the principal axis of rotation, and let OA, OB, be the 
two other principal axes at O. It is assumed that the top is maintained in constant 
angular rotation about the axis OC, which is regarded as executing vibrations of 
small amplitude about OZ, a fixed vertical axis. OA, OB, OC, are a system of 
rectangular axes which are principal axes fixed in the top; ZN is an arc drawn 
perpendicular to the quadrant BC produced and meeting it at N. The arc ZN is 
taken to be ® and the arc NC to be 9, both these arcs being supposed small. Let 
the constant angular velocity imposed on the top be 2; thus designating the three 
angular velocities of the axes OA, OB, OC, respectively as ,, ,, ®,, we have, 
since the arcs ®, 0, are by supposition, small and are rotating with constant angular 
velocity {2 about OZ at Z, 


o,= —(9 +2) along BC, . ; (1) 
o,=(P—O@) along CA, . ‘ (2) 
= in the direction AB. (3) 


G is the centre of gravity of the top, at the distance / from 0. 


3. DYNAMICS OF THE PROBLEM 


Let A, B, C, be the principal moments of inertia of the top at the point O. 
The inertia moments about the corresponding axes at O, fixed in the top, will be 


[Aw, +(C— B)o, 
or . (4) 
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about the, axis OA in the direction C to B; and similarly 
—[Bw,—(C — A), 
or — Bd—(C— A) 0?D+(A+B-C)Q0 . (5) 
about the axis OB in the direction C to A. Hence the equations of motion will be 
—ad —(c— b) 220 —(a+ b- sine,=0, (©) 
bb —(c— a) + (a+ b— 0) 20 + + mid QO? (7) 


where in engineers’ units W is the weight of the top, m=W/g, a=A/g, ¢tc., 
g being the acceleration due to gravity, and m/d®* the moment of dynamic unbalance 
(arising, for example, if OC is not precisely a principal axis), about an axis through O 
in the plane AOB. Let the axis of this unbalanced couple be inclined at an angle 

to OB in the direction A to B. The unbalance is assumed to be insufficient to 
have any appreciable effect on the principal moments of inertia of the top; it does 
however give rise to small products of inertia through the agency of which the 
unbalanced forces rotating with the impressed constant angular velocity become 
manifest as forced vibrations. 
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4. STEADY MOTION SOLUTION 
Clearly a particular solution of equations (6) and (7) presents itself of the form 


mld Q? sin <, 
[ (c—b)Q?- Wi] 
mld Q? cos 
which represents steady circular motion. 

The speeds for which the radius or amplitude of this steady motion become, 
in the absence of damping, infinite, are called whirling speeds. Thus if a, b, are 
both greater than c, then 0,, ®,, can never become infinite for any real values of ‘). 
There are therefore in such cases no real whirling speeds. If, however, {2 becomes 
infinitely large then 


mld sin <4 
tends to (c—b) (10) 


(c—a) 
If a, b, are both less than c, then there will be two real whirling speeds, namely 
2,2=WI/(c—b), 2,2=Wif(c-a) . . 


Clearly in this case if a>b then 2,2>,7. For the case a>c> b there will be 
only one real whirling speed, namely that given by 2,? above. 


5. STABILITY OF THE STEADY MOTION 


Next the stability of the steady motion is investigated, which as has been seen 
is conditioned by the lack of balance of the top. Ignoring this unbalance, a solution 
is assumed for the free vibrations of the form 


sin(wt+<*), P=, cos 
Inserting these values in the equations of motion, we obtain 
[aw? . (13) 
(a+ b—c) + [bw? . . (14) 
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On elimination of ©, and ®,, we derive the frequency equation, 


[aw? + (b—c) Q? + [bw* + (a— + WI (a+ (15) 
or 


abw' —{ [ab + (b—c)(a— c)] 2? — (a+ b) WI jw? + [(b - + (a- + WO =0 


For stability of the disturbed steady motion the two roots in w* of (16) must both 
be real and positive. Hence if there are two real whirling speeds, then apart from 
other considerations, the free vibrations are exponentially unstable for speeds 
between these whirling speeds; or if there is only one real whirling speed then 
stability is impossible for speeds above it. 


Furthermore { ab+(b—c)(a—c)} 2?-(a+b)Wl . ‘ (7) 
must be positive; and 


[ { ab+(b—c)(a—c)} 2?-(a+b) WI? 
— 4ab [ (b—c)(a—c)O* + (a+ b— 2c) + (18) 


must be positive; this latter condition may be expressed as 


[c(a+b—c)Q?— { (4ab/c)—(a+b)} WI]?> (19) 


6. THE SYMMETRICAL TOP 


For the ordinary symmetrical top b=a; in which case the condition for stability 
of the steady motion becomes 


Also if a> c there cannot be any real whirling speed; if however a<c there will 
be such a speed, given by 


and it may be shown that this will be above the minimum speed for exponential 
stability as given by (20). 


7. THE INVERTED TOP 


If the point of the peg is freely jointed, so ihat the top can spin suspended in an 
inverted position, then the preceding analysis is applicable to such cases provided 
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we change the sign of W. Thus we shall have for the steady motion controlle by 
the unbalance 


mld Q? sin ¢, 
(b—c)Q?- WI?’ 


22) 


mld Q cos 3 


If a, b, are both greater than c, there will be two real whirling speeds, naniely 


Clearly if in this case a> b then 2,?<,*. For the case a>c> b, there will be 
only one real whirling speed, namely that given by {2,? above. 


For the inverted symmetrical top the disturbed steady motion will be 
exponentially stable at all speeds; also there cannot be more than one real whirling 
speed and then only provided that a@ is greater than c. 
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THE POSSIBILITY OF THE DETERMINATION OF 


RATE OF CLIMB FROM ACCELERATION 
MEASUREMENTS IN LEVEL FLIGHT 
by 


E. C. PIKE, M.A., A.F.R.Ae.S. 


(Member of the Technical Staff of the Royal Aeronautical Society) 


SUMMARY 

Performance measurements for turbo-jet and turbo-prop aircraft are most 
conveniently reduced to standard conditions by the use of graphical methods. The 
case of rate of climb presents difficulties if readings are taken in the normal way 
during continuous climbs, because of the variation of weight, temperature, pressure 
and optimum climbing speed with height. 

Since rate of climb is proportional to the product of speed and acceleration 
in level flight, these quantities could be measured and used to determine the rate of 
climb at any height. 

The advantages of the method are listed and the technique required of the pilot 
examined. An automatic observer, including a longitudinal accelerometer, reading, 
in some cases, to within 0.001g, is necessary. 

The possibilities of error have been examined and have not been found to 
be serious. 

An example is given for a typical turbo-jet aircraft of the procedure in the 
determination of rate of climb at any forward speed for one value of N/ /@ and of 
W/s. To determine maximum rate of climb alone, fewer readings are sufficient. 
“ Non-dimensional ” curves of / 4 against N/ and W/« can be plotted directly 
from a series of such runs at different heights and different values of N. 


|. INTRODUCTION 
For the performance reduction of turbo-jet and turbo-prop aircraft, graphical 
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methods using the “non-dimensional” quantities V//0, N//0, W/o, 
v.//6@ etc. are the most convenient" 2), 


The level speed case is fairly simple because 


V//0=f,(N//0, for turbo-jets, 
) 


=9, (N//@, W/%, T,/9) for turbo-props, 


and during a level speed run at constant indicated height W/6, N//@ and 7,/6 
are approximately constant. Therefore for any fixed value of W/%, measured 
values of V//@ can be plotted against measured values of N//@ and 7,/6 (for 
turbo-props) or N/ 4/@ alone (for turbo-jets). The value of V/ ./@ under the required 
standard conditions can then be read directly from the faired curves through 
these points. 


The reduction of rate of climb results involves two difficulties, both of which are 
accentuated in the turbo-prop case. 


First, since rate of climb is a function of forward speed the additional variable 
V//6(=V,/4/@) is introduced into the equations, i.e. : — 


(2) 


J 0=f, (N/V 0, W/%, Vi/ /%) for turbo-jets, 
=, (N//0, W/o, T,/@, Vi//%) for turbo-props. 


In order to simplify the flying programme and the plotting of results, in practice 
the indicated air speed is varied with height in a pre-determined way so that the 
variable V;/ /@ is related to W/ and is effectively eliminated from the relations (2). 
With this procedure it is difficult to determine the variation of V, with height such 
that optimum rate of climb is achieved for all the values of N//6, W/ and T,/é 
which occur during the separate climbs. To determine the rate of climb at an 
arbitrary forward speed (not necessarily for maximum rate of climb) is even more 
difficult, and probably impossible. 


Second, in the usual method of measuring rate of climb, by observing times 
during a continuous climb, none of the variables in (2) remains constant. This 
means that preliminary cross-plots of the results are necessary before v,//@ can 
be plotted against the other “ non-dimensional ” parameters. However, it is possible 
effectively to measure rate of climb in level flight as the product of speed and 
acceleration, so eliminating the difficulties. 


* It has been assumed for simplicity of explanation that 74/6 is the additional variable in the 
engine performance of turbo-props (See Section 3 of Ref. 3). This assumption in no way 
prejudices the method. 
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The equation of motion for an aircraft can be written 


d 


(F—D)V= 


¥2+Wh) 


g dt + © 
= 
(F-D) = gi Gj). forclimb, hax 


dv . 


At any speed V the left hand side of equations (4) and (5) is the same whether 
the aircraft is climbing or not (assuming L=W, which is justified, particularly for 


turbo-jet or turbo-prop aircraft). Therefore if - &. is recorded in level flight, it 


can be equated to (1 which is equal to the rate of climb ft. / sec.) 


at the particular height (without taking into account the acceleration term ~ 4 c,) 


If the acceleration term is taken into account, as it should be for turbo-jet or 


turbo-prop aircraft, ‘ lA is readily calculated from runs at different heights. The 


value of V used in determining - a will usually be the optimum at each height 


(i.e. the speed for the greatest value of ve. or maximum rate of climb, at 


that height). 


The possibilities of using this method in practice are discussed here. 


Notation 
The suffix “a” denotes ambient or test conditions, the suffix 
standard conditions. 


“ 


s” denotes 


total aircraft drag (Ib.) 

propulsive thrust (Ib.) 

functions in equations (1) and (2) 
acceleration due to gravity (ft./sec.*) 
height (ft.) 
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LAS. indicated air speed (knots) 
L total lift (Ib.) 
N _ engine speed (r.p.m.) 
r.p.m. revolutions per minute 
T, absolute temperature in jet pipe (°K.) 
t time (sec.) 
V true air speed (ft. /sec.) 
Vax Maximum true air speed (ft./sec). 
V, equivalent air speed (knots) 
v. rate of climb (ft. / min.) 
W aircraft weight (Ib.) 
6 


ratio of absolute temperature to standard sea level value of absolute 
temperature 


® ratio of static pressure to standard sea level value of static pressure 
ao ratio of density to standard sea level density 
%,,¢_ functions in equations (1) and (2) 


2. ADVANTAGES OF THE METHOD 


2:1. For turbo-jets 


(i) Flying time can be saved by measuring rate of climb during the level speed 
run. If the aircraft is accelerated from about 0.5V nx. to 0.99V 4x. at a constant 
indicated height and constant N, measurements of rate of climb and maximum level 
speed can be made during the same run. 


(ii) During the run a “ partial climb” is measured at the same time, i.e. the 
rate of climb at any forward speed is known, and therefore the optimum value for 
the particular valués of N/./6 and W/® is ensured. This avoids making separate 
partial climbs to determine the optimum climbing speed as in the case of continuous 
climbs and therefore, the uncertainty is eliminated of the use of just one climbing 
speed at each height for the various values of N//@ and W/6 which occur during 
separate continuous climbs. 


In addition, the optimum climbing speed at each height may not give the best 
time to height, but results obtained by the suggested method can be used to determine 
immediately what is the optimum variation of forward speed with height to give 
the best time to height. 
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RATE OF CLIMB FROM ACCELERATION MEASUREMENTS 


(iii) During a run at constant N and constant indicated height, W/& and N// 6 
are constant. Therefore preliminary cross-plots of the results are not necessary, as 
they are for continuous climbs, where W/& and N//@ are varying continuously. 
This saves a considerable amount of time in the analysis and plotting of the results. 


(iv) By the use of the method suggested, allowance can be made for the 
acceleration term. 


(v) In obtaining v, from relations (4) and (5), the true rate of climb is determined 
directly. During continuous climbs the rate of change of pressure height is measured 
and has first to be multiplied by the factor 4,/, to obtain the true rate of climb. 


2.2. For turbo-props 


The advantages claimed above for the turbo-jet case apply with even greater 
force to turbo-props. This is because the appearance of the third variable T,/6 
in the engine performance equations complicates the flight testing and the analysis 
to an enormous extent. If continuous climbs are used, it is doubtful if the climbing 
speed variation selected is actually the optimum for all three variables and in 
addition a considerable amount of preliminary cross-plotting of the results is 
necessary, because of the continuous variation of all the “ non-dimensional ” 
parameters. 


3. PILOTING TECHNIQUE REQUIRED 


It is assumed that an automatic observer is fitted, recording N, 1.A.S., height, 
fuel contents, longitudinal acceleration and, in the case of the turbo-prop, jet pipe 
temperature. Outside air temperature is preferably recorded during the flight. 


The pilot is required to climb to the selected height, set the engine controls to 


‘the selected value of N (and to maintain the selected 7, for turbo-props) and then 


fly at constant altimeter reading during the acceleration from about 0.5V max. to 
0.99V max. The second and third runs (for different N) would be made at slightly 
different altimeter readings to allow for the fuel consumed during the flight. For 
instance, for the turbo-jet aircraft given in the example, the first run would be at 
20,000 ft., the second at 21,000 ft. and the third at 22,000 ft., so that W/@ is the 
same for each run. 


4 POSSIBILITIES OF ERROR 


The pilot may experience difficulty in maintaining constant indicated height 
during the inevitable changes of trim which will occur during the acceleration but 
it is not expected that the error would be serious. Other errors which may arise are: 

(i) Constant altimeter readings may not give a constant value of W/& due 
to fuel consumption during the run, and to position error. 
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HEIGHT MIN 


VARIATION 
(FT) sf 


Ss 06 09 1:0 


Fig. 1. 
Climbing effect. 


(ii) Increase of actual height at constant indicated height, due to position 
error, means that the aircraft is actually climbing very slightly. 


(ui) The longitudinal accelerometer may not be sufficiently sensitive. 


The significance of these errors is considered in what follows. 


4.1. Error due to variation of W/* during the run 


Two factors .affect W/& during the run; the reduction in W due to fuel con- 
sumption, and the change in & at constant altimeter reading due to position error. 
For aircraft with pitot heads having an altimeter correction which increases positively 
with speed, the effects of altimeter error and weight reduction act in opposite senses. 
Calculations for the typical turbo-jet aircraft of the example at 20,000 ft. and 
35,000 ft. show that between 50 per cent. and 90 per cent. of Vnax-, Which covers the 
climb region, a constant value of altimeter reading gives a value of W/® which is 
constant within 0.6 per cent. Since nearly all jet aircraft have leading edge pitot 
heads, for which the altimeter correction increases positively with speed, the 
assumption of constant altimeter reading is justified. 
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RATE OF CLIMB FROM ACCELERATION MEASUREMENTS 


4.2. Error due to shght climbing effect during run 


If constant indicated height is maintained during the run, and the altimeter 
correction increases positively with speed, the aircraft actually gains height. A 
slight error in the deduced rate of climb will result. Fig. 1 shows approximately 
the change of height for the typical turbo-jet aircraft used in the example. The 
error is seldom likely to be more than 6 per cent. and a correction can easily be 
made if necessary. 


4.3. Error due to sensitivity of the accelerometer 


As the example shows, a rate of climb of 3,000 ft. /min. at 220 knots corresponds 
to about 0.1g, so that for lower rates of climb the sensitivity of the accelerometer 
would have to be of the order of 0.001g. This accuracy has already been obtained 
with an improvised instrument, so that there should be no danger of error from this 
source, particularly as the development of accelerometers proceeds. 


5. EXAMPLE 


During a run at 20,000 ft. at 14,000 r.p.m. the following figures might be 
observed for a typical turbo-jet aircraft. The instrument readings have 
been corrected. 


N//@= cor =15,140 (from corrected value of r.p.m. and tempera- 


ture at time of test) 


W /a,= wap =54,400 (from fuel contents and altimeter reading) 
also o,=6,/6,=0.53. 
knots 158 178 497..: 247 237 257 276 287 329 


longitudinal 
acceleration: g 0.118 0.110 0.107 0.100 0.0928 0.0834 0.0652 0.0471 0 


‘ V; x 1.689 x 60 
and therefore, using v.=acceleration x — (0.53) 


v. : ft./min. 2590 2720 2940 3020 3060 2980 2500 1880 0 


These results are plotted in Fig. 2 and give a complete partial climb curve for 
N//@=15,140, W/S=54,400. To obtain only the maximum rate of climb it would 
not be necessary to plot as many points. ; 
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fewer test points ore required. 


iSO 200 250 
V; (KNOTS) 
Fig. 2. 
Example of results. 


For maximum rate of climb, in this case 


_ 3080 | 
v.//0= 0921 = 3340. 

The complete series of runs will give a “ non-dimensional ” set of curves from 
which the maximum rate of climb under standard conditions at any height and for 
any value of N can be obtained, and allowance can be made for acceleration during 
a continuous climb. 


A similar set of curves can be drawn giving the rate of climb at any arbitrary 
forward speed. 


The advantage of the fact that the “ non-dimensional” parameters are constant 
during the run is apparent during the plotting of the complete results for various 
heights and values of N. 
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A NOTE ON SUBSONIC AEROFOIL THEORY 
by 


JOHN W. MILES, Ph.D. 
(Department of Engineering, University of California, Los Angeles) 


SUMMARY 


A linearised theory for the subsonic, lifting surface problem is formulated in 
terms of Fourier integral solutions to Laplace’s equation. The symmetric and 
anti-symmetric problems of the first kind are solved explicitly, while the problems 
of the second kind depend on the solution of dual integral equations. The anti- 
symmetric problem of the second kind is cast in a variational form, from which 
certain well-known theorems may be deduced. 


1. INTRODUCTION 


The following development of linearised, subsonic aerofoil theory is derived 
from a recent analysis given for the supersonic problem“) which in turn was 
motivated by Schwinger’s variational approach to diffraction problems’: *). While 
the majority of the end results to be obtained are not new, their formulation, in 
certain respects, is rather novel, and the analysis is felt to be of interest per se. 


Earlier formulations of linearised, lifting surface theories are discussed by 
von Ka4rman“), and rather general approaches have been given more recently by 
Kussner“) and Reissner“). A Fourier transform approach of lifting line theory has 
been discussed by von Karman) and Fuchs"). 


Notation 
x co-ordinate of axis of rotation in pitch 
wing span 
reference value of chord 
drag coefficient 
lift coefficient 
rolling moment coefficient 


pitching moment coefficient 


Paper received April 1949. 
[The Aeronautical Quarterly, Vol. II, May 1956] 
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see equation (7) 

see equation (9) 

see equation (12) 

see equation (13) 

“the imaginary part of” 

unit vectors along x, y, z 
MacDonald’s function of order zero''?? 
MacDonald’s function of order one‘*?? 
perturbation pressure 

vector perturbation velocity 
integration operator. See equation (8) 
transform operator. See equation (4a) 
inverse of T 

x, y, Z components of perturbation velocity 


Fourier transform of u. See equation (4) 


P 
q 
R 
T 
T 
WwW 
U 
4 


free stream velocity 


x 
= 
N 


Cartesian co-ordinates (Fig. 1) 


RX 


wing incidence (z= —w in c) 

Dirac delta function; | (u)=1; (u)=0, 0. 

y co-ordinate of integration 

spectrum co-ordinates 

x co-ordinate of integration 

air mass density 

projeetion of aerofoil on z=0 

7 remainder of plane z=0 


@ velocity potential 


2. STATEMENT OF PROBLEM 


A thin wing is located near the xy(z=0) plane, its projection there being 
denoted by o (the remainder of the plane being denoted by 7). The theory is 
linearised in the sense that the boundary conditions are applied in o rather than at 
the aerofoil, and all perturbation velocities are assumed small compared to the free 
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}a(xy,o+) 


Fig. 1. 
Plan and section views of aerofoil. 


stream velocity, V, which is in the positive x direction. The perturbation velocity 
may be written 


q= (iu + jv +kw) V =VVo (1) 
V*o=0 
and the (linearised) perturbation (relative to free stream) pressure is given by 


where p is the mass density. It is required to find 9 (x, y, z) from a knowledge of 
the boundary conditions in the plane z=0 and at infinity. 


The boundary conditions require that the flow be tangent to the aerofoil and 
that the pressure be continuous except across the aerofoil. The top and bottom 
surfaces of the aerofoil will generally have different slopes, but the linearisation 
permits separate consideration of two problems consisting of the symmetrical and 
anti-symmetrical distributions of the incidence (measured in a section of constant y), 
corresponding to the thickness and camber (and twist) distributions, respectively. 
Hence, it suffices to consider only the half space z > 0, the boundary conditions 
being applied to u (x, y,0+) or w (x, y, 0+) (the case whose v (x, y, 0+) is specified 
is not generally of technical interest). 


In addition to the explicit boundary conditions in the plane z=0, it is required 
that the solutions be chosen such that ¢ will vanish far upstream, reduce to a 
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two-dimensional flow far downstream (the Trefftz plane)'*) and vanish as y, z —> 00, 
Moreover, the Kutta condition must be satisfied, ie. u(x, y, 0+) must vanish 
continuously at the trailing edge of c. 


3. SOLUTION 


Let lower and upper case letters represent functions of the co-ordinates 
(x, y) and their double, complex, Fourier transforms in the (u, v) spectrum, 
respectively, e.g. 


u(x, yy=T { U(u, v)} = du | dvU (u, vyexp [i(ux+vy)] . (4a) 
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U (u, { u(x, y)} = at | dnu(E, (4b) 


Then a solution to Laplace’s equation (2) which exhibits the desired behaviour at 
infinity (in the half space z >0) and. whose x derivative reduces to u(x, y) at 
z=0+ may be established as 


¢(x, y, { +78 exp [—(u? + v? Pz] U ¥) } (5) 


The branch cuts in the complex v plane are taken from +i|u| to +i, 
respectively, and the path of integration in the complex v plane may be taken 
anywhere in the strip | -4 (v)| <| provided z>0, as assumed. Proper 
behaviour for large values of y is guaranteed by giving -./(v) the same sign as y, 
with the real axis admissible for all y. The path in the » plane is taken along 
¥ (u)<0, and the contour may be closed in the upper or lower half of the 
u» plane for large positive or negative x, respectively, giving a vanishing solution far 
upstream and a potential solution (the contribution of the pole) far downstream. 
If the path is taken along the real axis in the » plane, the indentation under the pole 
#=0 yields equal and opposite contributions upstream and downstream, which are 
then matched by the contribution from the term 8 (u«)°) in such a way as to give 
the correct solutions both upstream and downstream. These choices determine 
essentially the arbitrary function of (x, y) introduced when u (x, y, z) is integrated 
with respect to x to determine ¢ (x, y, z). The Kutta condition does not enter the 
solution explicitly at this point since it is implicit in the specification of u (x, y). 


The upwash corresponding to equation (5) may be written, 
w (x, y)=9. (x, y,0+)=—T { G »)U (u, v) + (0+)] (6) 
G (u, (u? + 28 (u) | | ‘ (7) 
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The requirement z=0+ in equation (6) guarantees absolute convergence of the 
integrals in the limits ./ (u)= -/ (v)=0, but in the following the results will be 
formally interpreted for z=0 (cf. footnote ft). Applying the Faltung theorem‘? 
to equation (6) yields 


w (x, —(2)"' R{ g(x—, y—n)ué 
where R is an operator denoting double integration over all (£, ), and g(x, y) is 


the inverse transform of G(u, v). Carrying out the integrations required for the 


g(x, y)= —(iz | y| du K, («| y| )exp(iux)— y~* . (9a) 


= | du K, (u| y|)sin(ux)—y~? . (9b) 


where the path of integration for equation (9a) is indented under the pole »=0, 
and K, is a modified Bessel function of the second kind in the notation of Watson“??. 


The result obtained by substituting equation (9c) in equation (8) evidently has 
a direct interpretation in terms of a source distribution over z=0+, and, after 
integration by parts to give integrals possessing Cauchy principal valuest, it is in 
agreement with results obtained by rather more direct approaches. 


4. INVERSION OF INTEGRAL EQUATION 


Equations (6) and (8) represent a direct, closed form solution for w (x, y) if 
u(x, y) is known everywhere for z=0+. Conversely, if w(x, y) is known these 
results are integral equations of the first kind for U (u, v) and u(x, y), respectively. 
The integral equation may be inverted by taking the inverse transform{ of 


+ If z is kept finite (>0), it is found that 


g(x, y)=— {2 (2+y%)-" [1 +x 
and no difficulties occur at z=0, after integration by parts. 


tIn taking the reciprocal of G(, v), cf. equation (7), the delta function drops out since, in any 
event, G*(u, v) has a first order zero at »=0. See, for example, Ref. 10, p. 21. 
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equation (6) and applying the Faltung theorem, namely 
u(x, y= —T { G* (u, v) W (u, »)} 
u(x, y)= T { g* y—-n)wE, 
G* (u, Y= G~* (u? 


Taking the inverse transform of equation (12) yields 


g* (x, y)=(i/z) | . (13a) 


co 


= -(2/s)| uduK,(u|y|)sin(ux) (136) 
0 


Here again, the results are formal, convergence being obtained through an 
appropriate choice of paths in the transform planes and interpretation at z=0+, 
followed by integration by parts in (x, y). 


5. PROBLEMS OF FIRST AND SECOND KIND 


In a problem of the first kind either u or w is known over the entire plane 
z=0+, and the results of the last two sections therefore furnish exact solutions. 
In a problem of the second kind u and w are specified over different (but compatible) 
portions of the plane, and it is necessary to solve a “dual” integral equationt. 


Referring to the symmetric (s) and anti-symmetric (a) problems previously 
mentioned, w, and u, vanish over that portion of the plane (7) not occupied by the 
aerofoil («); hence, the determinations of either u, or w, from a knowledge of w, 
Or Ua, respectively, in « leads to a problem of the first kind, while the determination 
of u, or w, from w, or u,, respectively, in « leads to a problem of the second kind. 
The only exact solutions known for cases of interest in the latter category are the 
solutions for the circular"*®) and elliptic’) plan forms. 


6. UNIQUENESS 


While no adequate treatment of dual integral equations of the first kind appears 
to exist in the literature, they may be expected to exhibit characteristics similar to 


t See Ref. 11, pp. 334-342, for discussion of this type of integral equation. 
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other Fredholm integral equations arising in mathematical physics. (See, for example, 
Ref. 15.) In particular, eigensolutions may be expected for any given plan form, 
and these will generally possess singularities at the plan form edges. Various 
singularities are possible, but only singularities which are integrable, in the sense 
that they yield finite forces on the wing, can be admitted in the present problem, 
namely logarithmic and (—4) power singularities, in the symmetric and anti- 
symmetric problems, respectively"*). Now the particular solution to the integral 
equation may also be expected to exhibit a (—4) power singularity (in the anti- 
symmetric problem) at the plan form edges, so that the eigensolution may be chosen 
in such a way as just to cancel this singularity along the trailing edge, thus satisfying 
the Kutta condition and establishing the solution as unique.t The singularity along 
the leading edge remains and indeed is required (on the assumption of an ideal fluid) 
to eliminate drag (other than induced). 


The question of uniqueness is a recondite one, both mathematically and 
physically (no physical answer to the question can be given without considering the 
effects of viscosity in some detail), and it cannot be said that the question has been 
answered in any completely satisfactory way. It is, therefore, of some interest to 
remark that the same situation exists in other fields, e.g. diffraction’’), and has 
often been ignored. An interesting, if laconic, discussion of the supersonic case 
has been given recently by Hayes"*). 


7. VARIATIONAL FORMULATION 


A variational formulation of the dual integral equation of the first kind has been 
given by Schwinger”: *). Such a formulation of the present problem leads to a 
variational expression for the drag and is of interest only for the lifting surface 
problem of the second kind (angle of attack, i.e. downwash, distribution prescribed 
over aerofoil), since it is known, a priori, that the drag for a symmetrical aerofoil 
must vanish (ideal, irrotational flow). 


The drag is obtained by integrating the product of pressure, cf. equation (3), 
and downwash (w, (x, y)= — 2, (x, y)) over the aerofoil (c). The result ist 


R { [p (x, y,0—)—p (x, y,0+)] a (x, y)} 
=(4/c)R { (x, y) (x, y)} (14 


tit is also of interest to note that the discontinuity across the vortex sheet behind the wing is 
uniquely determined by the Kutta condition, or by any other uniqueness condition which 
might be formulated. 

{The matter of leading edge thrust (see, for example, Ref. 4, pp. 26-29) is passed over in silence 
here, except to state that the end results require no correction for such a thrust. It is not 
necessary to restrict the integration to «, since ua vanishes off o. 


29 


JOHN W. MILES 


It follows from Parseval’s theorem? that equation (14) may also be written 
Cyp=(4/o)R { v) A v)} 


where the bar denotes the complex conjugate. Alternative forms are given by 
substituting 2(x, y) or A(u, v) from the integral equations (6) or (8), respectively, 
to obtain 


p= R? { u(x, y—n)ulé,)} . (16) 


Cy=(4/7)R { U (u, »)G (u, v) U (u, } ‘ 


where R? indicates integration over both (x, y) and (€, »). If equations (16) and (17) 
are divided by the squares of equations (14) and (15), respectively, the results are 


Cy =(0/8x)[R { u(x, y)2 (x, y)} x { u(x, y)g y—n)ucé n)} (18) 
Cy-'=(¢/4)[R { U (uy, vA } xR (19) 


Equations (18) and (19) are variational in the sense that they are absolute 
minima (not merely stationary) with respect to arbitrary (not merely first order) 
variations of u or U about the exact solutions to the integral equations (6) and (8), 
respectively. A proof of this assertion, due to Schwinger, has been given 
elsewhere, it being specifically remarked that the present kernel is symmetric and 
positive definite. 


To prove the assertion that the kernel is positive definite, it may be observed 
that, since u (x, y) is real, its transform satisfies the relation 


whence U (un, v) U (n, v) is an even function of » and v, and only that part of 


G (u, which is éven in i.e. cf. equation (7), contributes to the result. 
In particular, equation (17) reduces to 


x 
(*)| vdv | U (0, v) |? 


0 


+Ref. 11, p. 51 (2.1.8). 
tSee Appendix IV of Ref. 19, the proof being for a vector integral equation with dyadic kernel. 
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as obtained (essentially) by von Karman‘), This result corresponds to the calculation 
of the induced drag in the Trefftz plane (x=00 or u=0). It may be remarked 
further that 


oo 


U (0, v)=(2z)-' dy exp (—ivn) | dé u(é, ») . 

and is therefore essentially the Fourier transform of the circulation about the 
“lifting line.” The result that the induced drag is independent of chordwise 


distribution of pressure (for fixed circulation) is of course the well-known Munk 
stagger theorem. (Ref. 4, p. 132.) 


The variational equivalent of equation (21) is deduced from equation (19) as 


oo 
U (0, v) 7 
(7), v)A (u, = 


The result that the reciprocal of Cy is a minimum, so that Cp is a maximum, with 
respect to variations about the true solution for the circulation is not new, for the 
usual (Fourier series) method of solving the Prandtl, lifting line, integral equation 
generally exhibits Cp in a positive definite, quadratic (in the expansion coefficients) 
form. See, for example, the work of Hildebrand’. 


8. LIFT AND MOMENT COEFFICIENTS 


To complete the specification of the forces acting on the wing, the lift and 
moment coefficients will be expressed in terms of the transforms. 


The lift coefficient is defined by 
R{ p(x, y, O—)—p(x, y, 0+)} 


Similarly, the pitching moment about x=a and the rolling moment about the line 
y=0 are referred to the coefficients 
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where c is a reference wing chord, and b is the wing span. These integrals may be 
expressed in terms of U (x, v), just as in the case of the circulation, cf. equation (22), 
since u, vanishes off the wing. The results are 


-(=)v@o 
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THEORY OF AN OSCILLATING SUPERSONIC 
AEROFOIL 


by 


GEOFFREY L. SEWELL, B.A. 
(King’s College, London) 


In the previous linearised theories of steady supersonic flow past a slender 
oscillating aerofoil, the shock line at the nose has been assumed to be a straight, 
stationary one, so that no vorticity is produced on passing through it. In the present 
treatment, the possibility that the shock may be both curved and moving is taken 
into account, but the result is exactly the same as that obtained by the previous 
treatments (e.g. Temple and Jahn"). 


Notation 


Uy, Pos Por Gov Mo constant values of flow velocity, pressure, density, speed of 
sound and Mach angle in the free stream 


Ox, Oy, Oz rectangular Cartesian. fixed right-handed axes. O mean 
position of nose, Ox direction of free stream, Oz direction of 
axis of rotation of aerofoil. 


time from a given fixed instant 


perturbation values of velocity, pressure and density 
respectively after passing through the shock 


transverse velocity of nose in direction Oy at time ¢ 
inclination of aerofoil to Ox at time ¢ 

angular velocity of aerofoil about Oz at time t(=9(t)) 
unit vector in direction Oz 

introduced and defined in context 

shock curve 

a point on § 

inclination of tangent at P to Ox at time f 

the line y=x tan u 


flow velocity components immediately behind P normal and 
tangential to the shock 


corresponding quantities immediately in front of P 
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e, (x, y, 0, d(x, 1) — introduced and defined in context 


y ratio of specific heat at constant pressure to specific heat at 
constant volume 


hy 


Shock 
Curve 


Fig. 1. 


Consider first the case of a flat aerofoil OA in a stream of constant speed u,, whose 
fixed direction will be taken to be that of the x-axis. Let Oy be the line perpendicular 
to Ox in the plane of oscillation, and let Oz complete a right-handed set of axes 
Ox, Oy, Oz. The aerofoil is oscillating with angular velocity » (t)(=6), while the 
nose, O, is vibrating along Oy with velocity v,(t). Let u,, p,, p, be the velocity, 
pressure and density respectively in the initial uniform stream, and let u,+u, 
Po +P. Po +p denote the corresponding gantities at (x, y, 0), a point in front of the 
shock. Moreover, since p, p, u, v,, , 6 are all small quantities whose squares we 
shall neglect, the equations of continuity, motion and energy become respectively 


(1) 


oP (SP +u, 22) 


where a, is the speed of sound in the initial uniform stream. Carrier) has shown 
that the solution of these equations may be written in the form 
u= grad H + curl (ck), 
p= Po Cr uy ox 


__ Po (OH =) 
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where k is the unit vector in the direction Oz, z=a(x—u,t, y) and A satisfies 


the equation 
2 2 ms 
(£ + Uy . : (4) 


The boundary conditions at the shock will now be determined. Let /,, m, be 
the normal and tangential components, respectively, of the velocity of flow relative 
to the shock immediately behind it at P; and let /, m represent the corresponding 
quantities immediately in front of the shock at P. The shock equations are 


m=m,, 
(e+ po) l=Pold, 
2y_ (P+Po 
y—-1\p+po 
y being the ratio of the specific heats of the gas. The solution of these equations is 
m=m,, 
2 
y+1 |, 


2 2 2 
p= y+1 Po a, ), 


I-1,= 


> 


Pit (l,,? —a,”) 


Since the shock is weak, /,—a, is small, and therefore, if we put /,=a,+« and 
neglect «*, it can be shown that 


m—m,=0, 


and consequently the shock condition becomes 
m—m,:1—1,:p:p=0: — ay: poy”: Py 


on S, the shock curve. The equation of S is y=xtanu,+d (x, t), where p, is the 
Mach angle for the initial uniform stream and d is small. It follows that, if 
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o(x, y, is also small, its value on S$ is equal to o(x, xtany,, f), correct to the 
first order of small quantities. This means that the value of # on S can be replaced 
by its value on L, the line y=xtany,. Thus, the shock condition is 


: on L . P (B) 


because Y—p, is a first order term, ¥ being the angle which the tangent to S at P 
makes with the x-axis. As p=a,*p on L, it is seen from (A) that f=0 on L. 
Therefore f(x—u,t, xtan u,)=0, and consequently f=0, since x—u,t and x tan 
are independent. The formulae for u, p and p may now be written 


u= grad H + curl (ck), | 


p=a,"p ar Uy Ox j 

The boundary condition at the aerofoil is that v=ox+u,4+v, on OA, and the same 
reasoning as that applied to the shock line proves that OA may be approximated 
tov=0. Thus 


v=oOx+U,9+v, on y=0. (D) 


The problem is now specified by (B, (C), (D) and (4), and it will be proved 
that Temple and Jahn’s solution satisfies all these conditions. Firstly, it satisfies 
(D), (C) and (4) because it is so constructed, and it therefore remains to be shown 
that it also satisfies (B). In the solution in question z=0 and H=0 on L. 


0H _ 


oH 
—tanu,=0 on L. 
y 


Therefore 
0} 


0 and re + 
Cx 


Thus “Rh <a te sin u, : cos, : 0 on L and therefore 
Cc 


Hence (B) is satisfied, which means that Temple and Jahn’s solution satisfies 
completely this linearised theory. 


Pp 4 


Fig. 2. 


On applying the above theory to a curved aerofoil, (D) is the only condition 
which undergoes any change. 
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Let the upper surface be of the form y=F (x) when the chord OA coincides with 
the x-axis. Thus PN =F (ON), and since the aerofoil is slender, F (x) is small. The 
inclination of the direction of the velocity at P to the x-axis is therefore 6+ F’ (x), 
so that the condition (D) is replaced by 


V=OX+U,0+, + UF’ (x) on y=0. 


Consequently, in order to obtain the complete solution for a curved aerofoil, we 
merely superpose the solution for a steady aerofoil of the same form on that for a 
flat aerofoil of the same length oscillating in the same manner. 
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THE LINEARISED THEORY OF CONICAL FIELDS 
IN SUPERSONIC FLOW, 


WITH APPLICATIONS TO PLANE AEROFOILS 
by 


S. GOLDSTEIN, F.R.S., F.R.Ae.S., 
and 
G. N. WARD 
(Department of Mathematics, The University, Manchester) 


SUMMARY 

In many important problems of supersonic flow, either for the whole field of 
flow or a part of it, the velocity components are constant on straight lines through 
a fixed point. Such velocity fields are called conical fields. In the usual theory 
of the linearised perturbations of a steady supersonic flow, the velocity is assumed 
to differ only slightly from a uniform undisturbed velocity, and in the defining 
equations and boundary conditions all non-linear terms in the components of the 
perturbation velocity (and their space derivatives) are neglected. In this paper the 
equations of linearised supersonic conical fields, and their general solution, are 
set out both for the region inside and for the region outside the Mach cone of the 
origin in the conical field. The results are applied to flow past plane triangular 
aerofoils with straight edges downstream of the vertex (of which there are six cases), 
to flow past those plane aerofoils of more extended shape and finite span for which 
the solution may be obtained by combining a finite number of conical fields, and to 
the problem of plane triangular vanes in semi-infinite free jets. In the applications, 
the velocity fields are determined in considerable detail, but a main purpose in 
setting them out was to exhibit the mathematical methods used and the physical 
considerations that enter in determining the mathematical solutions. 


Notation 
a,b,c,d __ in Sections 9, 10, 11, integers 


f,. fo. fy | arbitrary periodic functions of o +6 


arbitrary functions of ¢ 
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constant in Section 5 


modulus and complementary modulus of elliptic functions in 
Section 11 


modulus of elliptic functions, defined in equation (139) 
direction cosines of normal to surface of body 

in Section 11, defined by equation (148) 

in Sections 7, 9, 10, an integer 

in Section 7, an integer 


pressure of undisturbed stream; in equations (37), (38), (44) 
and (45), variable of integration 


perturbation pressure 

polar co-ordinates in (x,, y,) plane 

real variable defined in equation (21) 
complex variable defined in equation (60) 


values of ¢, being the limits of x, on aerofoils. (Equations 
(64), (110) ) 


perturbation velocity components 


4 
2 
3 


rectangular Cartesian co-ordinates. (Section 1) 


rectangular Cartesian co-ordinates after rotation of axes. 
(Section 3) 


x,,y, | rectangular co-ordinates in the plane z=cot «. (Equation (14)) 


A,, Az, A;, Ay, C, teal constants 


E(x),E incomplete and complete elliptic integrals of the second kind 
to modulus k 


complete elliptic integral of the second kind to modulus k, 


incomplete and complete elliptic integrals of the third kind, 
defined in equations (150) and (152) 


arbitrary periodic functions of « — 0 
arbitrary functions of t 
real constant in Section 6 


complete elliptic integrals of the first kind to moduli k, k’ in 
Section 11 


k, 
l, m,n 
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M_~ Mach number of undisturbed stream 
S., S», etc. | symbolic conical-field solutions in Section 12 
W velocity of undisturbed stream 
a angle of incidence, defined in Section | 
8., 8,, 8B, angles of sweep of leading or trailing edges of aerofoils 
arbitrary (small, positive) number, tending to zero 
complex variable defined in equation (23) 
polar angle in (x,, y,) plane 
variable of integration 
Mach angle 
angles defined in equations (13) and (47) 
density of undisturbed stream 
variable defined in equation (17) 


values of o at leading or trailing edges of aerofoils, defined in 
equation (40) 


perturbation velocity potential 


complex variable defined in equation (132) 
variable of integration 
“the real part of ” 


“ the imaginary part of ” 


l. INTRODUCTION. Linear perturbations of steady supersonic flow 

In the usual theory of the linearised perturbations of the steady supersonic 
flow of a gas, a velocity field is considered in which the velocity is supersonic 
everywhere. Viscosity, heat conduction, and heat radiation are neglected. The 
entropy of the gas per unit mass is assumed to be constant throughout the whole 
field. External forces, and therefore all variations in the potential energy per unit 
mass throughout the field, are neglected. The stagnation enthalpy (the sum of the 
enthalpy and kinetic energy) per unit mass, which is then constant along any stream- 
line, is assumed to be constant throughout the field, and the flow is irrotational 
(apart from the possible appearance of vortex sheets and singularities, as discussed 
in Section 6 et seq.). The velocity is assumed to differ only slightly from a uniform 
undisturbed velocity, and in the defining equations and boundary conditions all 


41 


i 
xX 
KH 
] 
|_| 


7 
4 
4 


S. GOLDSTEIN AND G. N. WARD 


terms are neglected that are of the second or higher degree in the components of 
the perturbation velocity and their space derivatives. 


Let (x, y, z) be rectangular Cartesian co-ordinates in the field, with the axis 
of z in the direction of the undisturbed velocity, whose magnitude is denoted by W. 
Let the undisturbed Mach number, equal to the ratio of W to the undisturbed 
acoustic velocity, be M, and let the Mach angle be u, so that u is the acute angle 
for which 


Denote the velocity components by u, v, W +w, and let » be the perturbation velocity 
potential, so that 


u=0¢/dx, v=0o/dy, w=00/ dz. (2) 


Also denote by p and p the undisturbed pressure and density of the gas, and by 
p+Ap the pressure at any point. 


From the equation of continuity, u, v, w, and ¢ all satisfy the equation 


Ox? + Gy? 


to a linear approximation and, from the integral of the equations of motion, 
Ap=—pWw, 


also to a linear approximation. 


We now assume that the disturbance in the flow is caused by the presence of 
a solid body, and, in specifying boundary conditions, we consider first the conditions 
at its surface. In some cases (e.g. for aerofoils) it is desirable to separate the effect 
of incidence, and therefore to begin with a precise definition of incidence. Generally 
the origin of co-ordinates will be taken as a point on the surface of the body such 
that no part of the surface is upstream of the origin, in z<0. (The origin is so 
chosen in all subsequent sections of this paper except at the end of Section 10.) For 
definiteness, the axes of x and y are assumed chosen so that they form a right-handed 
system with the axis of z, but may otherwise be chosen as most convenient. An 
incidence « is defined as a right-handed rotation about the x-axis through an angle 2 
from a specified standard position, and this definition is taken into account in 
choosing the axes of x and y; for example, if the surface of the solid body in its 
standard position is near a plane through the axis of z, this plane is taken to be 
y=0. In the standard position, let /, m, n be the direction cosines of the normal to 
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the surface drawn through a point of the surface into the fluid; then at incidence a 
the direction cosines are 


l, mcosa—nsinz, ncosz+msin 2. 


To accord with a linear approximation 2 must be small, and when only first powers 
of 2 are retained the condition that the normal velocity is zero at the surface is 


lu +(m—nz)v +(n+mz)(W +w)=0. (5) 


Since u, v, w are small compared with W, the multiplier of W must be small. The 
term in w and the term nav may be omitted, and the surface must be such that n 
is small. The boundary condition (5) becomes 


and, since 2 is small, it may be applied at the standard position of the surface, «=0. 


When the surface is near the plane y=0, and has small slopes (e.g. for thin 
aerofoils with flat or nearly flat boundaries), / and n are small, and 1 —m is small of 
the second order. The boundary condition (6) then simplifies to 


and may be applied on the part of the plane y=0 that is the projection of the surface. 


Other boundary conditions are derived from the consideration that with no 
part of the bounding surface upstream of the origin there is no disturbance for 
z<0, and u, v, w, and ¢ all vanish at z=0 (except perhaps at the origin and at any 
other points in the bounding surface in the plane z=0). Further, if no section of 
the bounding surface by a plane z=constant reaches to infinity, the disturbance will 
have reached only a finite distance for any finite value of z, and, in any plane 
z=constant, u, v, w, ¢, and all their derivatives, vanish beyond a certain finite 
distance from the axis of z. 


2. INTRODUCTION. Conical fields 


In many important problems of supersonic flow, either for the whole field of 
flow or a part of it, no fundamental length can be specified, and the velocity 
components are constant on straight lines through a fixed point. With this point 
as origin, u, v, w are functions of the angular variables only in a system of spherical 
polar co-ordinates; in Cartesian co-ordinates u, v, w are homogeneous functions of 
x, y, z of degree zero. Such velocity fields are called conical fields; their existence 
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was first recognised by Busemann"). When linear approximations are made, the 
mathematical analysis of such fields is particularly simple and elegant. There are 
other powerful and more general methods of considering solutions of the equations 
of Section 1; we are not here concerned with these other methods (to which some 
references are given at the end). The practical importance of conical fields, and the 
simplicity and elegance of their mathematics, entitle them to a separate detailed 
discussion. 


In the following sections we set out the equations of linearised supersonic conical 
fields and their general solution, and we apply the results to flow along plane 
triangular aerofoils with straight edges lying downstream of the vertex (which is 
taken as the origin), to flow past some plane aerofoils of a more extended shape, 
and to flow past plane triangular vanes in semi-infinite free jets. At zero incidence 
the aerofoils are in the plane y=0, so n=0 on both the upper and lower surfaces, 
and the boundary condition in (7) is simply 


on the part of the plane y=0 occupied by the aerofoil at zero incidence; the boundary 
condition is the same for both the upper and lower surfaces.* In such cases v is an 
even function, and u and w are odd functions, of y. 


Many results for plane aerofoils have been obtained by other authors, but 
the methods used here for obtaining both the fundamental relations for conical 
fields and the solutions for plane aerofoils are believed not to have been previously 
published in their entirety in the form here presented. The solutions for plane 
aerofoils seem to be more complete than those usually presented, and the results 
for plane vanes in a free jet are new. A list of references to the work of other 
authors is given at the end of this paper. Special mention should be made here of 
the work of W. D. Hayes, to which the attention of one of us (S.G.) was drawn when 
the main contents of this paper were presented in January 1948 in a lecture at the 
California Institute of Technology; Hayes’s methods have much in common with 
ours. The overlap with the works of other authors is much less. 


* The methods used in this paper may also be applied to more complicated examples, and a 
few such applications have been made by our co-workers. One such application is to the 
problem of the triangular aerofoil with finite dihedral and with edges swept back inside the 
Mach cone of the vertex. The aerofoil consists of two plane portions, which meet along 
the axis of z at an angle —28 and are each inclined, at zero incidence, at an angle 4 to the 
plane y=0. The boundary condition (6) is 

Fu sind+vcos 5=—Wacos on xsind+ycosé=0 for x 20 
for both the upper and lower surfaces. The conditions apply on the parts of the planes 
x sin8+ycos§=0 occupied by the aerofoil at zero incidence, and 5 has been taken as 
positive if the parts of the aerofoil are bent upwards (towards the axis of y) 
from the plane y=0. 
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3. FLOW OVER AN INFINITE PLANE YAWED WING WITH A 


STRAIGHT LEADING EDGE 


It will be convenient to set out here for reference the results for an infinite 
plane wing with a straight leading edge. 


When the leading edge is perpendicular to the undisturbed stream, it is taken 
as the axis of x ; the solution is independent of x, and u=0. The solution,* which 
is well known, is, for y > 0, 


v=w=0 for z<-—ycoth, v=—-Wa, w=Watany forz>ycoty, . (9) 
and, for y <0, 


v=w=0 for z<ycotn, v w=—Watanpy forz>-ycotpu. (10) 


Yz 
Fig. 1. 


If the aerofoil is yawed through an angle 6, (<47 — 4), so that its leading edge 
is along the axis of x’ (Fig. 1), the solution* is, for y > 0, 


u=v=w=0 for zcos 8, < ycot +x sin B,, 


v=—W2, w=Wacos 8, tan pn, 
for zcos 8, > y cot +x sin B,, 


* These solutions and the solutions in the following sections do not hold downstream of the 
trailing edges of the aerofoils. 
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and, for y <0, 
u=v=w=0 for zcos 8, <—ycot xsin 


v=—W2, w= —Wacos 8, tan», 
for zcos 8, >—ycot sin 


In this solution », is the acute angle defined by 
M cos 8, sin »»=1, cot? u.=M? cos? 8, — 1=cot? «cos? 8,—sin? . (13) 


and the incidence 2 is, by definition, a positive rotation about the axis of x, not about 
the axis of x’. The solution may be found by finding the relevant solution for 
depending on y and z’ only, and independent of x’, with the boundary condition 
v= —Waon y=0 for z’>0. Alternatively, we may note that the incidence 2 involves 
rotations «cos 8,, — «sin 8, about the axes of x’ and z’ respectively; the undisturbed 
velocity is compounded of W cos 8, along the axis of z’ and W sin 8, along the 
axis of x’. At incidence < the contribution to the component of the undisturbed 
velocity normal to the aerofoil is, therefore, Wz cos? 8, from the component W cos §,, 
and Wasin? 8, from the component W sin 8,, and their sum is Wz. The problem 
is therefore, to a linear approximation, the same as that of a wing with its leading 
edge along the axis of x’, placed in a stream with an undisturbed velocity W cos 8, 
normal to the leading edge, and set at an incidence (in the usual sense for a two- 
dimensional aerofoil of a rotation about the axis of x’) equal to z/cos £,. 


4. EQUATIONS AND GENERAL SOLUTIONS FOR LINEARISED 
CONICAL FIELDS 
Since, in conical fields, u, v, w are homogeneous functions of x, y, z of degree 
zero, if we write 


then u, v, w are functions of x, and y, only, and of » and @ only. 


In terms of r and @ the equations (of the form (3)) satisfied by u, v, and 
w become 


92 


For r> 1 we write the equation (15) for u as 


=| 
r (r?—1) r(r?>—1) 36?’ 
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and put 1/r=coso 


so tha do? * 


As r increases from | to 00, 7 increases from 0 to 47. 
The general solution for u in r> 1 is 
u=f,(¢+0)+F, —9), 
where f,, F, are arbitrary functions. Similarly 
v=f, +0)+F,(c —9), . (19d) 
w=f,(¢+0)+F, (c —9). (9c) 


Since u, v, w are single-valued, the fs and F’s must be periodic functions with 
period 2z. 


For r < 1 we write the equation (15) for u as 


ra-ry + - 


and put s= cosh~! (i [n= log 


O*u 
so that + 562 =(), 


Now let 


ptt 
Then for r <1, 

where .# denotes “ the real part of,” and g, (0) is an arbitrary function. Similarly 


When r=1, s=0 and e*=1; as r decreases from 1 to 0, 1/r increases from 1 to ©, 
s decreases from 0 to — 0, and e’ decreases from 1 to 0. Hence the interior of the 
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unit circle r=1 in the (r, 4) plane corresponds with the interior of the unit circle 
| | =1 in the ¢-plane; the boundaries correspond, and @ is the same in both planes, 


We have still to take account of the conditions that u, v, w shall be the 
components of the gradient of a function ¢, which itself satisfies equation (3). Since 
u, v, w are homogeneous functions of degree zero in x, y, and z, 

d (xu + yv+zw—9)=xdu + ydv+ zdw=0. . . - (26) 

If g, (Q=u+iu,, the component of the gradient of u, in any direction is equal 
to the component of the gradient of u in a direction at right angles, and similarly 
for v and w. Hence equation (26) implies that 

xdu, + ydv, +zdw,=0, 
and xg,’ O+ ye.’ 
Differentiation of (27) with respect to x, y, and z in turn leads to the equation 
[xg,” © + ye.” 


and two similar equations. It follows that 


, ot at 


By expressing derivatives with respect to x, y, and z in terms of derivatives with 
respect to r and @, it may be proved that 


a y-4(; | 
2z tan (s+i6)=(1—r*) C+ 
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and the general solution for u, v, and w inr <1 is 


u=-#{ 1g Qda, 


v= i| (A? - 1) g (A) da, 


w=2 tan | Ag (A) da, 


where g (A) is an arbitrary function of A. 


Since ¢ is a homogeneous function of degree zero in x, y, and z, 


from which it is easily checked that, if ¢=xu + yu + zw, 
then (u, v, w)=(0/dx, 2/2z)¢, and also (from (31)) that @ itself satisfies 


the equation (3). 


When r > 1, by use of the same arguments the relations 


xf,’ + 9) + yf.’ +9) +2f,’ +9)=0, 


xF,’ 0)+ yF,’ (o —9)+2F,’ —9)= 0, J 


are found for the derivatives of the functions in equation (19), and from these 
relations it is deduced that 


=cos (o + sin(o +9) : —tan p, 


=cos(o—6): —sin(o—6) : —tan 
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Hence the general solution for u, v, and w in r> 1 is 


6 


f (p) cos p dp+ { F (p) cos p dp, 


o+6 


f (p) sin p dp— | F (p) sin p dp, 


w=-tanu f(p)dp—tann | F(p)dp, 


where f (p) and F (p) are periodic functions of period 27 such that 
22 2n 2n 
{ | pdp=| F(p)sin p dp=0, 
0 


and similarly for F (p). 


5. THE VELOCITY COMPONENTS OUTSIDE THE MACH CONE 

OF THE ORIGIN* 

In an application to the flow over a plane triangular aerofoil the boundary 
condition v= — Wz is to be applied on the part of the plane y=0 occupied by the 
aerofoil at zero incidence, which corresponds to a certain portion of the line y, =0 
in the (x,, y,) plane. This portion may or may not lie, wholly or partly, outside 
the circle r=1; in this section we consider the cases when it does. 


The part of the (x,, y,) plane for which r>1 corresponds in the (c, 4) plane 
to a strip between and o=47. The axis y, =0 corresponds to (for x, > 0) 
and to 6=z (for x, <0), so the boundaries in the (c, 4) plane corresponding 
to those parts of a plane triangular aerofoil that are outside the Mach cone of its 
vertex lie along 9=0 and 6=7- (and are repeated at intervals of 27). They extend 
from «=0 (r=1) if the aerofoil passes through the Mach cone, but certainly do not 
reach to c=47(r=°0). We consider first the case in which the boundaries extend 
from o=0 to on 6=0 and to c=c, on (0, < 7, <4r, see 
Fig. 2a). Since 
(x/z) cot u.=x,=rcos 0=sec cos 8, d (39) 
this case is that of a triangular aerofoil whose edges are outside and on opposite 
sides of the Mach cone and are swept back at angles 8,, 8,, with 8, << 47-1», 
B, (Fig. 25), and 
sec 7,=cot 8, cot u, seco,=cot cot 


* The part downstream of any point P of the right circular cone with vertex at P and axis in 
the direction of the z-axis, and with semi-vertical angle equal to the Mach angle « will be 
called the Mach cone of P. The Mach cone of the origin corresponds with the circle 
r=1 in the (x,, y,) plane. 
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Now near 


since there is no disturbance beyond some finite value of r. Hence, from (195), 


fe (7 +6)+F, (c —6)=f,’ +6)+F,’ —6)= fa’ +9) F,’ (o — 4)=0, (42) 


$0 f.’ (¢ + —6)=0 


and (0 + 0)= —F, -0)=C 


near «=47, where C is a constant. Now and are constant 
along the lines 7 +@=constant and « —6=constant, respectively. (These lines are 
the characteristics of the differential equation (18).) If lines «+6=constant and 
¢—6=constant can be drawn from any point P to «=4 without intersecting the 
cuts from «=0 to along 6=0 and from «=0 to along then 
equations (43) still hold, and v=0, at P. Thus v, and similarly u and w, vanish 
everywhere for r > 1 except in the regions labelled 1, 2, 3, 4, in Fig. 2a. 


At the cuts the fs and F’s may be discontinuous. Regions | and 3 may be 
reached from o=42 along lines *—4=constant without a cut being crossed, so 
F,(¢—6) does not change as these regions are entered; but the regions cannot be 
entered from along lines + 4=constant if a cut is not crossed, so f, +4) 
may change discontinuously. In order that v should be equal to — Wa on the cuts, 
f,(o +) must, in fact, decrease discontinuously by Wz as « + 6 decreases through oc, 
when region 3 is entered, and as «+4 decreases through +,+7 when region | is 
entered. Similarly, as regions 2 and 4 are entered, f, (¢ +6) does not change, but 
F,(¢—6) must decrease discontinuously by Wa as «—@ decreases through o,, and 


through o,—7, respectively. Hence v=Wa throughout all four regions. To find 


ot+6 
u and w, we use (37). For example, on entering region 3, f (p) sin p dp 
decreases discontinuously by Wz as o + @ decreases through o,, so 


sin lim | f(p)dp=Wz. 
«99 
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The discontinuous decrease in u is 


+e 


o 


cos lim { f(p)dp=Wecotr, . . (45) 


% =e 


and that in w is similarly —Watanucoseco,. Hence u=—Wacote, and 
w=Watancoseco, throughout region 3. The values in the other regions are 
found similarly. The results are as follows : — 


Region u v w 
2 Wacotc, Watanpcoseco, } . (46) 
3 —Wacote, —~Wa Wz tan cosec 
4 Wacoto, — Watan cosec o, 


In these results v is an even function, and u and w are odd functions, of y, 
as they should be. ' 


Now consider the corresponding regions in the (x,, y,) plane and in physical 
space. Since rcos7=1, «+4=k corresponds to the line r cos (k—@)=1, which is 
the tangent to the circle r=1 at @=k. Similarly e—6=-—k corresponds to the same 
tangent. Along any such tangent r decreases from °° to | and then increases again 
to 00, and the complete tangent therefore corresponds to two characteristics in the 
(c, 6) plane lying between s=0 and o=47. Along one part of the tangent @ 
increases from r=°O to r=1; along the other part 6 decreases from r=°O to r=1. 
These correspond to the characteristics in the (c, 4) plane along which @ increases 
and decreases, respectively, as o« decreases from 4 to 0, and therefore to the 
characteristics o+0=k, o-—@=-—k, respectively. Thus the characteristics 
o+@=constant correspond to those parts of the tangents to r=1 on which @ 
increases as the circle is approached, and the characteristics « — @=constant to those 
parts of the tangents on which @ decreases as the circle is approached. 


Now at zero incidence the aerofoil in Fig. 2b corresponds, in the (x,, y,) 
plane, to the part of the axis y,=0 between x,=—seco,=-—cot8,cotu and 
x,=sec o,=cot 8, cot u, and the foregoing discussion shows that the regions 1, 2, 
3, 4 in the (c, 6) plane (Fig. 2a) correspond to the regions 1, 2, 3, 4 in the (x,,y,) 
plane (Fig. 2c), bounded by arcs of the circle r=1 (corresponding to «=0), by 
y,=0 (corresponding to 6=0 and 6=7), and by the tangents from the ends of the 
portion of y,=0 corresponding to the aerofoil (i.e. from x,=-—seco,, y,=0 and 
X,=seco,, y, =0, respectively). 
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Now let 1, “, be the acute angles for which 


M cos f, sin 4.=1, cos? 8,—1=cot* cos? 8, —sin? 
(47) 


and Mcos sin »,=1, cot? cos? 8, —1=cot? cos* B, — sin? 


(48) 


Then coto,=sin 8, tan coseco,=cot u cos tan Mo, 


cot o,=sin 8, tan 4,, cosec o,=cot » cos tan »,, 


and the values of u, v, w in (46) may be written in terms Of Bo, Bi, Mos i. For 
regions 2 and 3 the values are as follows : — 


Region u v w 
2 Wasin 8, tan », W2cos tan », 
3 —Wasin tan —W2a Wacos 8, tan 
In regions 1 and 4 the values of u and w are minus the values in regions 2 and 3, 


respectively, and v= —Wa. Outside the regions 1, 2, 3, 4, u, v, and w all vanish 
forr> 1. 


These are the same formulz as for infinite plane wings yawed through angles 
B., —8,, respectively (Section 3): outside the Mach cone of the vertex of the 
triangular wing, the effect of the vertex is not felt, as we should expect. Moreover, 
it is correct that the tangent boundaries to the regions 1, 2, 3, 4 correspond to the 
plane “ shock waves” or characteristic surfaces from the leading edges in physical 
space, whose equations (to a linear approximation) are (cf. equations (11) and (12)) 


zcos 8,=xsin8,+ycotm, (50) 
For example, the tangent boundary to region 3 is the line 
x, cOso,+y, sino,=1, ie. x, sin 8,+y, cot u.=cot cos B,, 
which corresponds to the plane 
x sin 8,+y cot u.=zcos B, 


in physical space, and similarly for the other boundaries. 


The relation between the “shock waves” from the leading edges in physical 
space, and the tangent boundaries in the (x,, y,) plane, may be seen in more detail 
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by another method. In the first place the “ shock waves” are the envelopes of the 
Mach cones of points on the leading edge at zero incidence. The Mach cone of the 
point (A cot 8, 0, A) on the leading edge in x >0, for example, is 


and the envelope of this family of cones is the pair of planes 
zcos B,=x sin 8, +y cot po. 


Now (x,, y,) are the values of (x, y) on the plane z=cot u, and any one of the above 
cones for which 0 <A < cot » meets this plane in the circle 


(x,-—Acot . . (54) 


For A=0 (the vertex of the triangular aerofoil) this is the unit circle, x,7+y,’=1; 
for A=cot » it is the point (cot » cot 8, 0). For any intermediate value of A it is 
the circle with centre at (A cot 8,, 0) and of radius 1—A tan »; its centre is on the 
line joining the origin to (cot » cot 8,, 0), and its radius is such that it touches the 
tangents to the unit circle from that point. The tangents therefore correspond to 
the envelopes of the cones in physical space. Moreover the insides of the cones 
correspond to the regions 1, 2, 3, 4 together with the inside of the unit circle—i.e. 
to the region in which u, v, w are not zero. (The region bounded by the continuations 
of the tangents from their points of contact to infinity, plus the unit circle, 
corresponds to the insides of the cones whose vertices are on the continuation of 
the leading edge into z <0; cf. the end of Section 10.) 


If one edge of the aerofoil is outside the Mach cone of the vertex and the other 
is inside, either regions 1 and 2 or regions 3 and 4 are absent, and the results for 
the remaining pair are unaltered, with w=v=w=0 elsewhere for r>1. If both 
edges are inside the Mach cone, u=v=w=0 everywhere in r> 1. 


6. THE VELOCITY COMPONENTS, OUTSIDE THE MACH CONE 
OF THE VERTEX, FOR A PLANE TRIANGULAR AEROFOIL 
ENTIRELY OUTSIDE THE MACH CONE 


If the whole of a plane triangular aerofoil with swept-back edges is outside 
the Mach cone of its vertex, then both edges are outside the Mach cone, but on the 
same side of it (Fig. 3a; 8, <8, < 47-1). 
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The aerofoil surface at zero incidence corresponds, in the (c, 4) plane, to the 
portion of for which <o <c,, with the same definitions of and as 
before (equation (40)). Except in the regions 1, 2, 3, 4 in Fig. 3b, in which the 
dashed lines are parts of e+ 0=0,, c+0=0,, u, v, and w are zero as before. The 
regions 1 and 2 are separated by 6=0, 0 < « <c,, and on the boundary the normal 
velocity and the pressure must be continuous, i.e. v and w must be continuous. 


As before, when the region 3 is entered, f, (¢ +) decreases discontinuously by 
Wa as «+6 decreases through o,, and F,(¢—) does not change. When the 
region 4 is entered, F,(«—6) decreases discontinuously by Wa as «—@ decreases 
through o,, and f,(¢ +6) does not change. Thus v is constant and equal to —Wa 
throughout regions 3 and 4. wu and w are found as before, and their values are the 
same as in regions 3 and 4 in (46). 


In passing into region 1 from region 3, F, (c —@) does not change; let f, (¢ + 4) 
increase by K as «+6 decreases through 7,, so v becomes —-Wa+K _ throughout 
region 1. By the argument as before, u will increase by Kcotc, and w by 
~K tan »coseco,. Hence in region 1 


u= —Wacote,+Kcotc,, v= —Wa+K, w=tan (Wacoseco,—K coseca,). (55) 
In passing into region 2 from region 4, f, ( +) does not change, and F, (¢ — @) must 
increase by K in order that v should be continuous on the boundary between regions 
land 2. In region 2, therefore, 


u=Wacotec,—Kcote,, v=—Wa+K, w=-—tanp(Wacoseco,—K (56) 


Since. w must be continuous on the boundary between regions 1 and 2, w=0 in 
both these regions and 


With the same definitions of », and », as before (equations (47) ), equations (48) 


still hold, and in terms of 8., 8,, “, #, the values of u, v, and w in region 1 are 
as follows : — 


— Bo) ta jo Wa(1- w=0. . (58) 


cos tan cos 


The values in region 3 are exactly as in (49). In regions 2 and 4, the values of 
u and w are minus the values in regions 1 and 3, respectively (w=0 in region 2), 
and the values of v are the same as those in | and 3, respectively. 


The regions 1, 2, 3, 4 in the (x,, y,) plane are shown in Fig. 3c . 
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There is a discontinuity in u across the boundary between region 1 and 
region 2, so there is, downstream from the trailing edge, a vortex sheet of uniform 
strength (which is continued inside the Mach cone; see Section 9). 


7. THE VELOCITY COMPONENTS INSIDE THE MACH CONE OF THE 
VERTEX 


For r< 1, u, v, and w are the real parts of functions g, (0, g,(Q, and g, (0) of ¢ 
(equation (24)). From the results of Sections 5 and 6, u, v, and w are known on the 
circle r=1 in the (x,, y,) plane, and therefore on the circle | § |=1 in the (-plane. 
The part of the aerofoil, at zero incidence, inside the Mach cone of the vertex 
corresponds to the whole or part of the axis y,=0 between x, = +1, and therefore 
to the whole or part of the real axis in the (-plane inside the circle | (|=1; along 
this part of the real axis in the (-plane, v= —Wa. Hence, for |¢|<1, the real 
part of g, (0 is known on the unit circle |¢|=1 and on any cut there may be 
inside it, and we may find g, (() at all points inside the unit circle in the (-plane. 
If g, (©) has no singularities for |¢| <1 and ¢ not on the cut,* the solution for 
g, () is unique. There are, however, cases in which the solution that is free from 
singularities inside the unit circle off the cut is not physically acceptable, and we 
must return to the question of the occurrence of such singularities; we may say at 
once that their occurrence is connected with the appearance of trailing vortex sheet, 
(Sections 9, 10, 11.) When g, (¢) has been determined, v is known; the relations (31), 
together with the values of u and w at any one point on the unit circle, enable u 
and w to be determined (and their values at other points on the unit circle may 
then be checked). 


To solve for g, (0) it is advantageous to transform from the (-plane to a f:plane 
by the conformal transformation 


cos (1—r’} sin 0 


From (23) ; = 


(60) 


The circle r=1 and its interior correspond with the circle |¢|=1 and its 
interior; the circles are transformed, in the f-plane, into the upper and lower sides 
of the parts of the real axis for which | ¢ | > 1 (Fig. 4); their interiors are represented 
on the whole of the f-plane with cuts along the portions | t | > 1 of the real axis, 
the upper halves (0 < @ < =) being represented on the upper half-plane ( 7 (t) > 0) 


* The singularities on the cut are at its ends and satisfy the conditions listed, in particular 
condition (ii). 
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and the lower halves (—~7<6<0) on the lower half-plane (./ (t)<0) 
(.£ denotes “the imaginary part of.”) The axis x,=0 and the imaginary axis in 
the (-plane, inside the unit circle in-each case, correspond with one another, and are 
transformed into the imaginary axis in the f-plane. 


(at i oO) 


¥,) -PLANE 


re) 


YO 
(at—i 
Fig. 4. 


The point (x,, 0) on the x,-axis inside the unit circle in the (x,, y,) plane 
corresponds (from (60) ) with the point t=.x, on the real axis in the t-plane. For a 
point on the unit circle in the (x,, y,) or (-plane, f=sec 0; at 9= +0, t=sec o + 0i, and 
sec > is the abscissa, in the (x,, y,) plane, of the intersection with the x,-axis of the 
tangent to the unit circle at 0= +o. Thus the point f=.x, on the real axis in the 
t-plane corresponds with the point (x,, 0) if | x,|<1, and with the points of 
contact of the tangents to the unit circle from (x,, 0) if | x, | > 1. 


Now v= — Wa on the part of the x,-axis on which the aerofoil at zero incidence 
is represented, i.e. on this part of the x,-axis, v= — Wa at (x,, 0); and, if | x,|>1, 
it follows from Section 5 that v= — Wa also at the points of contact of the tangents 
from (x,, 0) to the unit circle. Hence v= — Wa on the real axis in the f-plane for 
all points with abscisse between the extreme values of x, (i.e. of x/(z tan »)) on the 
surface of the aerofoil at zero incidence. 


Since u, v, w are the real parts of functions of ¢, they are also the real parts 
of functions of t. We put 


u= RG, (t), v= w= AG, (1) 
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and from (31) and (59) it follows that 


where (1—7*) is the branch, in the cut ¢-plane, that is positive when f¢ is real 
and < l. 


Essentially the process we shall follow is to solve first for G,’(¢). Certain 
restrictions may be immediately placed on G,’ (t) in order that the solution may 
correspond to a physically possible flow :— 


(i) On |§|=1, v is zero at (= +i and in a neighbourhood of each of these 
points; will be regular near (= +i, and these points correspond with 
t=0o, Since 


G,’ (t) must be O(t-*) at most as t—> oO in order that g,’(¢) should be finite 
at (= +i. 


(ii) In order that the lift force on the aerofoil should be finite, the pressure, 
and therefore (from (4))w, must be integrable over the aerofoil surface, the values 
being equal in magnitude and opposite in sign on the two sides of the aerofoil. 
Over the part of the aerofoil outside the Mach cone of its vertex, w is finite. Over 
the part inside the Mach cone w, if not finite, must at least be integrable. It follows 
from this argument, and from the correspondence of the line x/(z tan u)=x, with 
the point t=x,, that w (= ZG, (t)) must be integrable along the real axis in the 
t-plane between the extreme limits of x, on the aerofoil. In the sequel we shall see 
that G,’(f), and therefore also G,’(t), may have a factor (t—1,)'+* or (t—1,)", 
where f, is one of the limits of x, on the aerofoil, and m and m are integers. The 
preceding argument shows that m must be >-— 2, and m must be >- 1. 


(iii) In all cases u, v, and w must themselves be finite at x,= +1. We shall 
see that G,’(t) may have factors (1 +7)'+*, where n is an integer; G,’ (#) and G,’ () 
then have factors (1+7)'~*, and nm must be > 0. 


We shall consider the various cases in some detail, giving formule for the 
functions G, (t), G, G, in the upper half-plane (y, > 0) in each case. The 
velocity components u, v, w, for r<1, are the real parts of these functions; in only 
one case do we write out explicitly expressions for these real parts in terms of x, 
and y,. We shall not mention again the values in the lower half-plane (y, < 0), 
since they are covered by the fact that v is an even function, and u and w are odd 
functions, of y,. The values of v on the part of y,=0, | x,|<1, corresponding 
with all or part of the aerofoil surface at zero incidence, and its values on r=1, are 
used in obtaining G, (¢), and it may be verified from the final formule for G, that 
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they have been correctly used. It may also be verified that u and w are continuous 
on r=1, the values found in the following sections agreeing, on r=1, with those 
found in Sections 5 and 6. For this purpose it will be convenient to record the 
values in Sections 5 and 6 in terms of the notation of the following sections. In all 
the following sections we shall have 


t,=seco,=cot 8, cot », t,=seco,=cot cot . (64) 
Then, for ¢, > 1, 


t, (t.2— 1)~* tan u=tan cosec o,=cos tan . (66) 


(see equation (48) ), and similarly with suffix 1. The formule in (49) for the velocity 
components in regions 2 and 3 of Fig. 2c (for which ¢, > 1, t, > 1) become 


Region u v w 
2 Wa(t,?— Wa, Wat, 1)~* tan (67) 
3 —Wa (t,?- 1)-', = Wa, W at, (t,? — 1)-* tan 


The formule in (58) for the velocity components in region | of Fig. 3c become 


The result most often required is the pressure on the aerofoil surface. From 
equation (4), Ap= —pWw, and to a linear approximation the value may be taken 
on y=0 instead of on the surface of the aerofoil at incidence. For each case we 
shall give the value of w on the relevant part of the upper side of the plane y=0, 
designated as y= +0. The value on the lower side is equal in magnitude and 
Opposite in sign. 


It will appear that, in all the cases in which, at zero incidence, the aerofoil passes 
through the z-axis, G,’ (t) has no singularity at =0, so (from (62) ) G,’ (0)=0, and 
it is a general result that on either side of the aerofoil w, and therefore also the 
pressure on the surface, reaches its extreme value at x,=0 in such cases. 


8. BOTH EDGES OUTSIDE THE MACH CONE OF THE VERTEX, ON 
OPPOSITE SIDES 
For this case (Fig. 2), for y,>0, u, v, and w vanish for r>1 except in 
regions 2 and 3, where their values are given by (67). Hence, on the real axis in 
the t-plane, 
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(69) 


v= AG, ()= —Wa for -t, <t<t, 


for and t<—tf, 


The solution for G, (¢) may be written down immediately. In the upper half-plane 


to 


G, (t)=- * tog 


(70) 


with arg (¢—t,)=arg (t+1t,)=0 for ¢ real and > 
Then 


Wa(t,+t,) 
ix (t—t,)(t+t,) 


G,’ 
(so the conditions (i), (ii), (iii), in Section 7 are satisfied), and from (62) 


We tan 


Hence (with G, =G,=0 at t=00) 


1 l l+t,t } 


on the upper half-plane. 


By taking real parts we find that in the upper half-plane (y, > 0) and for r< L, 


{ sin 8, tan cos ( _ sin 8, tan cos~' (76) 


Wa _, 1+¢,x, 
{ cos 8, tan Cos ) A, tan p, } (78) 
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=[(t.-x,? 
B 


wf and .&. are positive, and the inverse cosines are between 0 and =. 


The value of w on y= +0, | x, | = 1 is given in (49) and (67) for the regions 
where it is not zero. For | x, | <1, on y= +0, 


1- 1+t,x, \. 


9, BOTH EDGES OUTSIDE THE MACH CONE OF THE VERTEX, AND 
ON THE SAME SIDE 


For this case (Fig. 3), for y, >0, u, v, and w vanish for r>1 except in 
regions | and 2. Their values in region | are given in (68), and in region 3 in (67). 
On the real axis in the f-plane 


v=0 for t<—1 and fort>t, 
=—Wa for iot>t, 


1) 
t, (t.7—1) 


Hence G,’ (t)=0 on the real axis for | ¢| > 1,¢ or 


= 


The pressure, and therefore also w, are odd functions of y, and are continuous 
across the axis y,=0 inside r=1. Hence w=0 on y,=0, | x,| <1, ie. w=0 on 
the real axis in the f-plane for |t|<1.* It follows that Z G,’ (t)}=0, and (from 
(62)) that G,’ (t)=0, on the real axis in the t-plane for |¢| <1. Hence, on 
the real axis, 


RG,’ (t)=0 for |t| >1, t At, or t,, FG,’ for |t|<1,. (82) 
and G, has finite discontinuities at t, and t,. We therefore write 


Git), . . (83) 


G,' ()= 


where a and b are integers and G (f) is a rational function of t, real when f is real. 


* The component of the gradient of v normal to the real axis is also zero. 
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The residues at the poles ¢, and ¢,, multiplied by zi, must be equal to the discon. 
tinuities in G, at those points along the real axis. 


If G()=O (t") as t—> then from condition (i) in Section 7, m> a+ b +1; 
from condition (iii), a0, b>0. Hence m>1; G has a singularity, which is at 
least a pole; if it is a pole, a=b=0, m=1, and G is of the form constant /(1-t,) 


We found in Section 6 that the boundary between regions 1 and 2 contains a 
vortex sheet of constant strength. The vortex lines coincide with the streamlines, 
and so, to the present approximation, may be taken as straight lines parallel to the 
axis of z, and the strength of the sheet is then, with the same approximation, 
independent of z. It follows that the whole sheet is of constant strength, and is 
bounded by the streamline through the vertex of the aerofoil, i.e., to the present 
approximation, by the axis of z. Hence u must have a finite discontinuity at r=0, 
so G must have a pole at r=0 and be equal to A,/t, where A, is a constant. (That 
the singularity in G must be at t=0, and cannot be a higher singularity than a poke, 
also follows from the consideration that the pressure, and therefore w, must be 
continuous along the real axis in the ¢-plane except at ¢, and ¢,.) The value of A, 
is found by equating the residue of G,’ (t) at t=t, to iWz/z. Hence 


G,’ (t)= to (t.—t,) (1 


(t,2—1) 


Wa (t.—1,) 1 
ix 


From (62), G,’ ()=- 


Wa t, (t.—t,) tan pu 


Hence, on the upper half-plane, 


G,(=- in t.)— log (t—1,) +“ . (87) 


t,-t t, (t.7—1) 


(88) 


G, (j= Flog . ‘ ; (89) 


The vortex sheet is of constant strength 


2 (t,—t,) 
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There is a logarithmic infinity in v at ¢=0, i.e. along the approximate edge of 
the vortex sheet, the axis of z. The linear approximation is therefore not adequate 
in a neighbourhood of that axis. 


10. ONE EDGE OUTSIDE AND ONE EDGE INSIDE THE MACH CONE 


Let the edges be swept back at angles 8,, where 8, <P, 
<4iz+yu. We have to distinguish two cases: (i) 8, > 42 (Fig. 5); (ii) 8, <4 
(Fig. 6). In the first case the aerofoil extends through the z-axis; in the second case 
itdoes not. When 8, >47 the flow in the vicinity of the edge swept back at the 
angle 8, resembles that near a subsonic leading edge, and singularities may be 
expected there on a linear theory. If, however, 8, <4, the flow near the edge £, 
is like that near a subsonic trailing edge, and we apply the Kutta-Joukowski 
condition that the velocity shall be finite there; there will be a trailing vortex sheet 
of uniform strength downstream from this edge, with its other boundary on the 
z-axis. In both cases the velocities in the region 3 are those in (49) or (67). 
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The boundary conditions in the ft-plane are the same in both cases; the 
differences in the solutions come from differences in the singularities. By the same 
arguments as in Section 9, on the real axis in the f-plane 


v=—Wa fort, <t<t, | 
=0 for t<—1 and for t>t, - 


w=0 for -l1<t<t,. 


RG, (t}=0 for t<—1 and for t>t,, 


FG,’ (t)=0 for -l<t<t, 
on the real axis, and G,’ (t) is of the form 


G, =i . (94) 


where a and b are integers and G (f) is a rational function of f, real when ¢ is real. 
From the second and third of the conditions in Section 7, b >~—2 and a> 0; from 
the first condition, if G (t)=O (t-") as t—> ©, m>a+b+2. 


In the first case, when 8, > 47, t, is negative, the wing passes through the 
z-axis, and there is no trailing vortex sheet. G (t) will have no singularities, and 
m=0. Hence a=0, b=-—2, m=0, and G is a constant A,. The value of A, 
is determined by the discontinuity in v at t=t,; the residue of G,’(t) at its pole 
t=t, must be iWa/x. Hence 


_ iWa(t,-t,p? (1+0 


(95) 


Wa 1 


[2 


Hence, on the upper half-plane, 
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(100) 


On y, = +0, for t; <x, < 1 (ie. for the part of the aerofoil surface, at zero incidence, 
inside the Mach cone) 


_2Wa —cot B, (t,—t,)(1—x,) 1) (x, 


cot 8, and ft, being negative. 


In the second case, when 8 < 4x, 1, is positive, and there is a trailing vortex 
sheet. G,’ (t) is still given by equation (94), but the Kutta-Joukowski condition, that 
the velocity shall be finite at t=/,, requires that b=>>—1. Hence m>1, and G (#) 
has a singularity. By the same argument as in Section 9, the only singularity in G 
can be a pole at t=0, and a=0, b= — 1, G (t)=A,/t. The constant A, is determined 
in the same way as A,, and 


Wa to (t.—t,) 1 


G,’ 7 (t, +1) t(t- 


iW2t(t.-th 
+1 


(t, +1) (t—t.)(t—t,¥ 


Hence, on the upper half-plane, 


(105) 


(106) 
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On y,= +0, for tp < x, <1, 


(t.— 1). - 4) 


2W2 
8, tan tan=? [ 


The vortex sheet is of constant strength 


y 


It is convenient at this stage to exhibit a solution for the tip of a swept-back 
aerofoil. Let O’ (Fig. 7a) be the tip of an aerofoil whose vertex is O, and let 00’ 


z 


Fig. 7a. Fig. 7b. 


be a leading edge swept back at an angle 8,<(4x-—y. With O’ as origin, this 
edge appears as a supersonic leading edge swept forward at an angle 8,. Let the 
other edge through O’ be swept back at an angle 8,, where 4x—n < 8, <4r+u. 
(The cases in which 8, and 8, >4x+, may be dealt with similarly by 
an extension of ‘the analysis of Sections 8 and 9.) Thus one edge through 0’ is 
outside the Mach cone of O’, and the other edge is inside it. We are for the present 
concerned to find the velocity field inside the portion of the Mach cone of O’ that 
is not intersected by the Mach cone of O. 


The velocity over the aerofoil between the two Mach cones is given by the 
formulz in Section 5, and, in the upper half-plane, by the formula for the region 3 
in (49) and (67). In the (x,, y,) plane, with O’ as origin, these formule give u, 2, 
and w in the region 3 of Fig. 7b, and in particular over its curved boundary. Hence, 
with ¢, defined as before (equation (64) ), the expressions for u,v, and w for region 3 
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in (67) give the values on the upper side of the real axis in the f-plane for t >t, 
and ¢ <—1 (since the relevant portion of the unit circle in the (x,, y,) plane trans- 
forms into these parts of the real axis in the f-plane). Over the part of the unit 
circle in the (x,, y,) plane that is not the boundary of the regions 3 and 4 the velocity 
components vanish, so, on the real axis in the f-plane, w=v=w=0 for l1<t<t. 
On the part of the aerofoil itself (taken at zero incidence, as usual) inside the Mach 
cone, v= — Wa, so if 


t, =cot 8, cot », (110) 


v=—Wa on the real axis in the t-plane for —1<1<t,. Furthermore, for 
t, <t< 1, w=0 on the real axis, by the same reasoning as in Section 9. 


There are again two cases to be considered. If 8, < 47, t, >0, the edge 
swept back at an angle £, is a subsonic leading edge, at which we expect singularities 
on a linear theory, and there is no vortex sheet. (This is the case illustrated in 
Fig. 7.) If, however, 8, > 42, t, <0, then the edge is a subsonic trailing edge, 
at which we apply the Kutta-Joukowski condition that the velocity shall be finite, 
and there is a vortex sheet extending (according to the linear approximation) along 

=0 from x, =f, to x,=0. 


The boundary conditions for G,’ (t) on the real axis in the t-plane are now 
RG, (t)=0 for t<t, and fort -7G,'()=0 fort,<t<1. (ill) 


The solution is found in the same way as before. For the first case (8, << 47, t, > 0, 
no vortex sheet), on the upper half-plane, 


(112) 


tant 


(114) 


On y,= +0, -1 <x, (/e. on the part of the upper aerofoil surface inside 
the Mach cone of O’), 


t, 


_ 2Wa {cot B, 


\i+t, +c0s 8, tan tan” 


(115) 
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For the second case (8, > 47, t, <0, vortex sheet from x,=f, to x,=0), on 
the upper half-plane, 


1 
} 


_ _ 2iWat, tan 


On y= +0, for -1 <x, 


w= —— 
us 


(t.+ I(t, = 


-1 
cos Ae tan 


The strength of the vortex sheet is 


11. BOTH EDGES INSIDE THE MACH CONE 


Let the edges be swept back at angles 8, and f, (Figs. 8a, 8b), where 
<8, 4x and 8, <8, All the velocity components, u, v, and w, 


in 
2Wa (t. —t, (1—t)]! 

a 
t! 
t 
z 2 | 
Fig. 8a. Fig. 8b. 
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are zero on the circle r=1 in the (x,, y,) plane, and therefore also on the real axis 
in the ft-plane for |t|=1. Also v=—Wa on the real axis in the f¢-plane for 
t,<t<t., where f., t, are defined in (64), and, by the same reasoning as in 
Section 9, w=0 for -1<t<t, and for t, <t<1. Hence 


(121) 


(t)=0 for t<-1, for t, t <t,, and for t> 1, 


(t)=0 for -1 <t <t, and for t,<t< 1, 


on the real axis, and 


where a, b, c, d are integers and G (f) is a rational function of f, real when ¢ is real. 


There are again two cases. In the first (Fig. 8a), 8, > 47, t, <0, the aerofoil 
at zero incidence extends through the z-axis, both edges are subsonic leading edges, 
at which singularities are to be expected on a linear theory, and there is no vortex 
sheet. In the second (Fig. 8b), 8,<{47, t, > 0, the aerofoil does not extend 
through the z-axis, the edge swept back at an angle £, is a subsonic trailing edge, 
at which the Kutta-Joukowski condition (that the velocity should be finite) is applied, 
and there is a trailing vortex sheet extending (according to the linear approximation) 
along y,=0 from x,=?f, to x,=0, corresponding with the real axis in the ¢-plane 
from to 


In the first case, G (t) has no singularities; in the second it has a pole at r=0. 
The conditions on G,’ (t) in Section 7 lead to the conclusions that in the first case 


a=b=0, c=d=-2, G()=A, .- (123) 
and in the second case (since u, v, w must be finite at t=1,), 
a=b=0, c=-2, d=-1, G(W=C,/t, 
where A, and C, are real constants. 
Hence in the first case 
A,(i-?} 


(t- 
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and, from (62), 


iA, 


GO= 


iA,t tan p 
whereas in the second case 


C, (1 


(128) 


(129) 


iC, tan u 


(130) 


The constant A, and C, are found, after integration, from the condition that 
in each case, on the real axis, v=0 for |¢| > 1 and v= —Wa for t, <t<t,. 


We now integrate for the first case (8, > 47, t, <0, no vortex sheet). G, (¢) 
and G, (¢) are expressible as elementary functions, but G, (¢) is an elliptic integral. 
Let 


_ 


2 (t.—t,) 


= 


(131) 


k, k’ being positive, and let K and E be the complete elliptic integrals of the first 
and second kind, respectively, to modulus k, and K’ the complete elliptic integral 
of the first kind to modulus k’. We define x by means of Jacobian elliptic functions: 


(1 —t,) 
(1—t.) (t—1,) 


(t,—t,) (1 
° 


dn (x, ° (132) 


and write E (x) for the elliptic integral of the second kind, 


sn (x, k)= cn (x, k)= 


E()= | dn? (6, k)dé. 
0 
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dx 1 (1—t,)(1+t.) 


dE(x)_1f +) 
aa 


The values t=1, f,, t;, — 1 correspond with x=K, 0, iK’, K +iK’, respectively. 
On the upper half-plane the integral for G, that vanishes at t=1 is 


(136) 


When t, > ¢ > ¢,, the last term on the right, x, and E (x) are all purely imaginary, so 


v= — Waa 37) 


~ 20 +4) + @ 


and A,= 


It may now be checked that v=0 when 1 is real and | ¢|> 1. 


The expression 2E—(k’)?K may be transformed by a transformation due to 
Gauss. A new modulus &, is introduced, where 


and the variable of integration y in the defining integrals for E and K, 
0/2 2 


E= | sin’ yp dy, K= | sin’ dy, 


0 


is changed by writing 


_ (1+k)siny ,_ sin’ y,) 


and it may be proved that 
2E—(k YK=(1+k)E,, 


where E, is the complete elliptic integral of the second kind to modulus k,. 
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Integrals for G, and G, on the upper half-plane whose real parts vanish when 
t is real and | t| > 1 are 


G,()= + t- }. 


Hence, on y, = +0, for t; <x, < ¢, (¢, being negative), 


In the special case when the aerofoil is symmetrical about the z-axis, t, = —1., 


and 


k?=1-t,3, 
1+f, 


In this case the expression (145) becomes 


W ct,” tan u 
(147) 

We turn now to the second case (8, < 42, t, >0, a vortex sheet on y,=0 
between x,=/, and x,=0), and integrate equations (128), (129), and (130). Again 
G, (¢) and G, (¢) are expressible as elementary functions, whereas G, (¢) is an elliptic 
integral. To integrate G,’(t), k, K, E, x, and E(x) are introduced as before 
(equations (131), (132), (133), (140)). In addition, let 


_ 


m t, (1+t.)’ (149) 


so that 


and 1 >k?>m>0. Also introduce the elliptic integral of the third kind 


dy 
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for which 


dF 1 (1 +1.) ¢—1,) 
x} (151) 


This elliptic integral has not been tabulated (Fletcher, Miller, and Rosenhead, An 
Index of Mathematical Tables, London, 1946, p. 318), but is expressible in terms of 
theta functions. The corresponding complete elliptic integral will be denoted by F, 
so that 


K 


0 


do 


Then on the upper half-plane the seieaad for G, that vanishes at t=1 is 


G, (p= { -1)[K-F-x+F } - 


_ 1-?f 


(t.—t,) 


When t, >t >t,, the last term on the right, x, E(x), and F (x) are all purely 
imaginary, so 


to—t, (1+t)(1—2,) { E+(k?/m-1)(K-F)}° 


and (155) 


It may now be checked that v=0 when ¢ is real and | t| > 1. 


Integrals for G, and G, in the upper half-plane whose reai parts vanish when 
tis real and | ¢| > 1 are 


Hence, on y, = +0, for t, <x, <t., 


2C, tan (x,—t,\* 


w= 
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The strength of the vortex sheet is 


2x |C, | 


As 8,—8,—> 0, t,—-t, —> 0, m—> 1, k—> 1, E—> 1, (k?/m—1)(K —F)—-0, 


_ 2C, tan Wat,tanp Wacos 8, 


and 


(160) 


Also _ 
t.—x, 


(161) 


where l,=zcos B.—xsinB., . (162) 


As 8,.—8,—> 0, sin 8./sin 8, —> 1, and the aerofoil, at large distances from the 
vertex, of the order of (8,—,)~*, becomes an aerofoil of infinite span yawed at an 
angle 8, > 47—,; on y=0, /,/I, is the ratio of the distance of a point from the 
trailing edge to its distance from the leading edge; and the formula (158) for w on 
y, = +0 for t, < x, <t, becomes the usual formula for an infinite yawed plane wing 
of constant chord. The strength of the vortex sheet tends to zero with f,—f,. 


12. FLOW OVER MORE GENERAL PLANE AEROFOILS BY FINITE 
COMBINATIONS OF CONICAL FIELDS 


The solutions in Sections 8, 9, and 10°may be combined to give the flow over 
certain extended classes of plane aerofoils at small incidences. We shall consider 
here only aerofoils the flow over which can be obtained by the combination of a 
finite number of conical fields. All supersonic leading edges, and all subsonic edges, 
of such aerofoils must be straight. The occurrence of a straight subsonic edge sets 
a limit to the extent of such an aerofoil; such an edge is inside the Mach cone of 
its upstream end, and a conical field (or finite combination of conical fields) cannot 
be found with the downstream end as origin to give correctly the flow over an 
adjacent edge on the downstream side. More extended classes of aerofoils, including 
aerofoils with curved leading edges, which may be obtained by superposing con- 
tinuous distributions of conical fields—i.e., by integrating conical-field solutions—are 
not considered here. 


The method of combination may be explained by considering the plane aerofoil 
shown in Fig. 9. Fig. 9a shows the aerofoil at zero incidence, in the plane y=0; 
to a linear approximation it is also the projection on to the plane y=0 of the surface 
of the aerofoil at a small incidence z. OP, OQ are straight supersonic leading edges: 
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PR, QS are straight subsonic edges; OR, OS, which meet PR and QS in R and §, 
are the traces on the plane y=0 of the Mach cone of O; PU, QT are parts of the 
traces of the Mach cones of P and Q, respectively; and RTV, SUV are parts of the 
traces of the Mach cones of R and S, respectively. 


As before, since u and w are even functions, and v is an odd function, of y, 
we need consider the solution only for y=>0. The regions a, b, c, d, e, f, etc., are 
all regions in y>0. The region a is the part of the Mach cone of O that is outside 
the Mach cones of P and Q; the region b is bounded by the Mach cones of O and P 
and the plane characteristic surface through OP, and the region c is defined similarly 
with Q in place of P. The region d is the part of the Mach cone of P outside the 
Mach cone of O, and similarly for e with Q in place of P. The region f is inside 
the Mach cones of O and P, and outside the Mach cones of R and Q; g is defined 
similarly with Q, S, and P in place of P, R, and Q, respectively. The region h is 
inside the Mach cones of P and Q, and outside the Mach cones of R and S. Sections 
through the regions are shown in Figs. 9b-9f. 


The velocity fields in the regions a, b, c were found in Sections 5 and 8 (for 
the effects of the breaks in the leading edges at P and Q are confined to the interiors 
of the Mach cones of these points). The fields in the regions d and e are conical 
fields with origins at P and Q, respectively, of the type calculated at the end of 
Section 10 (since the effect of the break in the leading edge at O is confined to the 


inside of the Mach cone of O). The fields in the regions a, b, c, d, e are therefore 
found as in previous sections. 


Now denote by S,, S», S., Si, . . . the solutions for the velocity components in 
the regions a, b, c, d,. . . , respectively, with the following stipulations. In S, and 
S., the values of the velocity components are constant, so §, and S. may be considered 
to be defined at all points. The solution §, may be continued throughout the whole 
of the region in which it would apply for an infinite triangular aerofoil, which is the 
interior of the Mach cone of O, and may be considered defined at all points of that 
region. Similarly S, and S. may be continued throughout the whole of the Mach 
cones of P and Q, respectively, and may be considered to be defined in those regions. 
(It is to be understood, however, that S,, S;, . . . , though now defined in more 
extended regions, are solutions for the given aerofoil only in the regions a, b,. . . , 
respectively). We shall now call the fields deduced from S,—S, and S.—S, the 
“complementary conical fields” with P and Q, respectively, as origins, and shall 
show that the solution in a more extended region inside the Mach cone of O may 
be obtained by subtracting these complementary fields from the field S,. In fact, 
the solutions for the regions f, g, h are given by 


S; =S, Sa), 
S, =§, (S. or: S.), 
Si =§, (Sp Si) (S. —S.). ) 
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These solutions clearly satisfy the approximate (linearised) field equations. They 
also satisfy the boundary condition v= — Wa on the aerofoil surface at y=0, and 
they make the velocity components continuous over the common boundaries of f 
with a and d, of g with a and e, and of h with f and g. The value of w, and therefore 
of Ap, on y= +0 in the regions f, g, and A may now be immediately deduced from 
the values for the solutions S,, S,, S., Sua, Se. 


The Mach cones of R and S include portions of the plane y=0 that are not on 
the projection of the aerofoil surface, and on which, to preserve continuity of pressure, 
w=0. This condition cannot be satisfied by combinations of the above type that 
also satisfy the other boundary conditions, and so no such combinations are valid 
anywhere in the Mach cones of R and S. Extended regions of the aerofoil surface 
may be treated only by superposing continuous distributions of conical fields—i.e., 
by integrating conical-field solutions. 


The foregoing solution yields the pressure over the whole aerofoil surface if the 
trailing edge of the aerofoil is any line (not necessarily straight, or in straight 
segments) that lies in the region bounded by OP RV SQ O in Fig. 9a and whose 
angle of sweep nowhere exceeds 47 —u. 


The solutions found in this and the preceding sections may be used to calculate 
the resultant aerodynamic forces and moments on a given aerofoil. The lift force, 
parallel to the y-axis, and the moment about the z-axis, may be calculated by 
integrating the pressures over the upper and lower surfaces. At a subsonic leading 
edge, however, the velocity on a linear theory becomes infinite; if r, is normal distance 
from the edge in the aerofoil plane, the velocity ~ Dr,~* as r, —> 0 (where D is a 
constant); there is a suction force on the leading edge, acting normally outwards from 
the edge and in the aerofoil plane, whose magnitude, per unit length of the edge, is 


(1 —M? 


. (164) 


where 8 is the angle of sweep of the edge. (See Refs. 8 and 27.) This suction must 
be included in the calculation of the drag and side forces and yawing moment on the 
aerofoil. 


13. PLANE VANES IN FREE JETS 

In this section we solve the problem of a plane triangular vane in a semi-infinite 
free jet, which is bounded by a plane vortex sheet through the vertex of the: vane 
and normal to the plane of the vane at zero incidence. The gas outside the jet is 
assumed to be at rest relative to the vane, with pressure equal to the free-stream 
pressure in the jet. 


With the origin of rectangular Cartesian co-ordinates at the vertex of the vane, 
the jet is taken to be in the half-space x > 0, with its boundary on x=0. The vane 
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is at a small incidence z; when «=0, the vane is in the plane y=0. The undisturbed 
gas velocity in the jet is W along the axis of z. 


On a linear theory the displacement of the jet boundary, as well as the incidence 
of the vane, may be neglected in formulating the boundary conditions. The pressure 
must be continuous on the jet boundary, so, on a linear theory, w must be zero there, 
The boundary conditions are that on y=0, x > 0, 


(165) 


v= — Wz on the projection of the vane, 


and w=0 elsewhere, 
and also w=0 on x=0.. ‘ (166) 


We shall consider only vanes having one swept-back leading edge outside the 
jet, and the other inside the jet, so that the vane passes through the vortex sheet, 
x=0. We shall solve for the disturbance velocity (u, v, w) in the half-space x > 0. 
In the (x,, y,) plane, since v is an even function, and u and w are odd functions, of 
y,, it will suffice to find the solution in the upper right-hand quadrant. 


Let the leading edge inside the jet be swept back at an angle §,, and put 
t,=cot 8, cot» as usual. We have two cases to consider: one when 8, < 47~1, 
t, > 1, and the leading edge is outside the Mach cone of the origin, and the other 
when 42 > 8, > 4x>—y, 0<t, <1, and the edge is inside the Mach cone. 


The (x,, y,) plane is shown (for x, > 0) for the first case in Fig. 10. For rr >1, 
the velocity components are zero outside region 3, and in region 3 they have the 
same values as in (67). We note that the condition w=0 on x,=0 is satisfied 


for r>1. 
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To find u, v, and w for r< 1, we transform to the f-plane, and find the solution 
in the upper right-hand quadrant of that plane. Consider the boundary conditions 
in the ‘-plane. From the values on r=1, we have that u, v, and w vanish on the 
real axis for ¢ > ¢t, and have the values in (67) (region 3) fort, >¢> 1. From (165), 
y= Hence, in particular, 


(167) 


for 


=(0 for t>t, 


on the real axis. The condition w=0 on x,=0 becomes, for r<1, AG, (t)=0 
on the imaginary axis. Hence, on the imaginary axis, -7 G,’ (t)=0, and, from (62), 
F G',(t)=0, so v is constant; since v=0 when t=i00, v=0 everywhere on the 
imaginary axis. 


The boundary conditions for G,’ (t) are 


(t)=0 on the real axis for t>O(t 
(t)=0 on the real axis for (t+ * (168) 


(t}=0 on the imaginary axis. 


The first condition in Section 7 still applies, and G,’ (t) is O (t-*) at most as t—> ©. 
The discontinuity at t=f, in the value of v along the real axis shows that G,’ (¢) 
has a pole there with residue iWa/x. There is also a pole at the origin, since 
v= —Wa2 when the origin is approached along the real axis and wv is zero on the 
imaginary axis; the residue at the origin is —2iW2/z. G,'(t) has no other 
singularities in the upper right-hand quadrant of the t-plane or on its boundary. 


The solution for G,’(t) and G, (ft) with these conditions is unique. G,’ (ft) 
is an odd function of t, the solution being 


2 


W 2_ 42 


From (62) 


2Wat,” 
at (t? 


G,’ (t)= 


2W 2t,” tan 
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whence, in the upper right-hand quadrant, 


G, ()=- { (173) 


It may be verified that the boundary conditions on u, v, and w are all satisfied. 


On y= +0, for0< x, <1, 


4 


t.?—x, 


This solution may be used in combination with other conical-field solutions to 
give the flow over certain vanes of finite span, as explained in Section 12. 


Both in this solution, and in the solution following for the second case, there is a 
logarithmic infinity in u at t=0, corresponding with the axis of z, where the vane 
crosses the boundary of the jet. The linear approximation is therefore not adequate 
in a neighbourhood of that line. 


In the second case, when 42 > 8, > O0<t,< 1, v=w=0 along the 
imaginary axis in the f-plane as before. Along the real axis, the known values are 
now 

v= for 0<t<t, 


=( for t> 1, 
u=w=0 ior t> 
Thus, on the real axis, 4G.’ (t)=0 for 0<t and fort > 1, 
JF (t)=0 for t.<t<1, 
and on the imaginary axis, .7 G,’ (t)=0. J 


In the upper right-hand quadrant the solution for G,’ (f) (satisfying the relevant 
conditions in Section 7) and for G, (f) is 


2Wet,> 
it 


G.()= (7-5) —tanh-! |}. (179) 
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2Wat,* 
nt (t.? 


G,’ (t)= 


2W tan 


{1-4 tanh-* (1- }, 


G,’ (t)= 


2W at, tan u t 
(t,2— 


G; (= 


and on y= +0, for0 <x, <h, 


2Wcot B, 
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